Applications of a tight-binding total energy method for transition and noble metals: 
Elastic Constants, Vacancies, and Surfaces of Monatomic Metals 
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A recent tight-binding scheme provides a method for extending the results of first principles cal- 
culations to regimes involving 10^ — 10^ atoms in a unit cell. The method uses an analytic set of 
two-center, non- orthogonal tight-binding parameters, on-site terms which change with the local en- 
vironment, and no pair potential. The free parameters in this method are chosen to simultaneously 
fit band structures and total energies from a set of first-principles calculations for monatomic fee and 
bcc crystals. To check the accuracy of this method we evaluate structural energy differences, elastic 
constants, vacancy formation energies, and surface energies, comparing to first-principles calcula- 
tions and experiment. In most cases there is good agreement between this theory and experiment. 
We present a detailed account of the method, a complete set of tight-binding parameters, and results 
for twenty-nine of the alkaline earth, transition and noble metals. 



I. INTRODUCTION 

Recently, a total energy tight-binding (TB) method 
was introduceda wherein the parameters were fit to the 
total energies and band structures obtained from a lim- 
ited set of first-principles calculations. The method pro- 
duces very good structural energy differences, elastic con- 
stants, phonon frequencies, vacancy formation energies, 
and surface energies for the transition and noble met- 
als. In this paper we thoroughly examine the method for 
these elements, and for the alkaline earth metals. 

Section |l| describes the method, including the ana- 
lytic expressions for the tight-binding parameters and the 
scheme for fitting them to first-principles theory. The 
remainder of the paper discusses the predictions made 
by this method. Sec. Ill shows that the method cor- 



rectly predicts the ordering of several metallic phases, 
including the 29 atom aMn phase, and correctly pre- 
dicts the energy differences between the lower phases. 
Sec. [V discusses the calculation of elastic constants. 

-|Can also be performed by first- 
we compare our TB results to 



Since these calculatic 
principles techniques ,B 
those obtained from first-principles to assess the accuracy 
of the method. We can also use the method to study 
systems which have relatively large (20-200 atom) unit 
cells, a regime which requires laborious first-principles 
calculations, but which can be done very quickly by our 
TB method. In Sec. ^ we discuss the structure and en- 
ergetics of vacancies, and in Sec. Vl we consider surface 
energetics. Finally, in Sec. VII we summarize our results. 



discuss possible methods to improve our agreement with 
first-principles theory or experiment, and consider the 
course of future improvements to the method. 



II. THE TIGHT-BINDING METHOD 

Density functional theory (DFT),J using the Kohn- 
Sham ansatz for the kinetic energy,Q tells us that the 



total energy of a system of electrons moving in a solid 
can be written in the form 



i?[n(r)] =^/(Ai-e,)e. + i^Kr)] 



(1) 



where n{r) is the electronic density, Si is the Kohn-Sham 
eigenvalue of the i*^ electronic state, /i is the chemical 
potential, and the sum is over all electronic states of the 
system. Since we will be primarily interested in studying 
metals, we take the function /(/i — e) to have the Fermi 
function formfl 



/(^) 



(2) 



where /3 — l/{kT). Typically we take T between 2 and 
5 mRy. The functional _F[rt(r)] contains the remaining 
part of the DFT total energy: the ion-ion interaction en- 
ergy, the parts of the Hartree and Exchange-Correlation 
energy not included in the eigenvalue sums, and correc- 
tions for double counting in the eigenvalue sums. 

Many tight-binding methods use (0) to provide a nat- 
ural separation between the parts of i?[n(r)] which can 
be treated by tight-binding, and the parts which must be 
treated by other means. El t3 In these cases the eigenvalue 
sum is calculated by tight-binding methods, and F[n(r)] 
is approximated by a sum of pair potentials. 

We have developed, an alternative method of applying 
tight-binding to (|^),a based on the fact that the Kohn- 
Sham method allows an arbitrary shift in the potential. 
If this shift is defined to be 



where 



Vo = F[n{r)]/Ne , 
7Ve = ^/(M-£.) 



(3) 
(4) 



is the number of electrons in the system, then the eigen- 
values Si are each shifted by an amount Vo, to the new 
values 



1 



e[ = e, + 14) . 
The total energy (|l|) then becomes 



(5) 



(6) 



where /i' = /i + Vb is the shifted chenwcal potential. 

By the density functional theoremja the shifted eigen- 
values e[ can be considered to be functions of the crys- 
tal structure, including volume, primitive lattice vectors, 
and internal parameters. A tight-binding method which 
reproduces the e[ over a range of structures will then 
solve the total energy problem (||) or (H) without resort 
to an additional term. !_. 

The two-center Slater-Koster formulatiorO of tight- 
binding with a non-orthogonal basis breaks the problem 
into the calculation of three types of parameters: on-site 
parameters, which represent the energy required to place 
an electron in a specific orbital, Hamiltonian parameters, 
which represent the matrix elements for electrons hop- 
ping from one site to another, and overlap parameters, 
which describe the mixing between the non-orthogonal 
orbitals on neighbor sites. The eigenvalues £\ can be de- 
termined once the parameters are evaluated for a given 
structure. The basic method we use is to give the Slater- 
Koster parameters simple algebraic forms, with param- 
eters chosen to reproduce first-principles results over a 
wide range of structures. 

In an orthogonal tight-binding calculation, applying 
the shift (H) to each eigenvalue would be equivalent to 
shifting each diagonal element of the Hamiltonian ma- 
trix by an amount Vq. In the non-orthogonal case the 
effect is slightly more complicated. It is clear, however, 
that the Hamiltonian and overlap parameters should not 
directly depend on the shift Vq. Thus the effect of (||) 
can only be accounted for by changing the on-site pa- 
rameters. These parameters must now be sensitive to the 
local environment around each atom. We describe this 
local environment by introducing a "density" associated 
with each atom, defined by 



= ^exp[-Ali?,,]^;(i?,,) , 



(7) 



where i (j) denotes the type of atom on site i (j), A is 
a parameter which will depend on the atom types, Rij 
is the distance between atoms i and j, and Fc{R) is a 
smooth cutoff function. 



i^e(i?) = {l + exp[(i?-i?o)/£]}- 



(8) 



which we use to limit the range of the parameters. In the 
calculations shown here we typically take Rq — 14.0 Bohr 
and £ ~ 0.5 Bohr, which effectively zeros all interactions 
for neighbors more than 16.5 Bohr apart. Depending on 
the structure and lattice constant, this radius will include 
from 80 to 300 neighboring atoms. 



Although in principle the on-site terms should have 
off-diagonal elements due to the overlap of the ophite 
wave functions with neighboring atomic potentials ,113 we 
follow traditional practice and only include the diagonal 
terms, keeping only the terms corresponding to s, p, and 
d orbitals. The on-site term for atom i is 



2/3 



CiaPi + diaP^ , (9) 

where a = s,p, or c?, and pi is given by (0). To show the 
correspondence between this theory and previous theo- 
ries, note that p in (Q) has a form similar to a pair poten- 
tial. Thus if (|) were linear in p, we restricted ourselves to 
an orthogonal Hamiltonian, and kept hia independent of 
the angular momentum, then (|^) would just be a method 
of parametrizing _F[n(r)] by a simple pair potential. Seen 
in this light, the density dependent parts of on-site term 
@ can be regarded as a generalized pair potential. 

In general the Hamiltonian and overlap Slater-Koster 
parameters may depend upon the structure and local en- 
vironment of the atoms, even in the case of the two- 
center approximation. This approach, while valid, will 
lead to many complications when we extend the method 
to include other crystal structures and binary alloys. We 
therefore constrain the form of these parameters so that 
they depend only on the distance between the two atoms. 
Both the Hamiltonian and overlap parameters are as- 
sumed to have the same functional form. 



P^(i?) = (e^ + f^R)eM-9lR]Fc{R) 



(10) 



where 7 indicates the type of interaction (e.g. ssa, pdn, 
etc.), R is the distance between the atoms, and Fc{R) is 
given by (||). 
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FIG. 1. Values of the two center Slater-Koster integrals for 
the Hamiltonian matrix as a function of distance for molybde- 
num, obtained from the tight-binding parameters ( |lo| ) 
integrals are labeled according to the standard notation! 

In the monatomic case, restricting ourselves to s, p, 
and d orbitals, there are three on-site Slater-Koster pa- 
rameters (^), with thirteen parameters required to com- 
pletely specify them. Similarly, there are ten Hamil- 
tonian and ten overlap parameters (p^, each requir- 
ing three parameters. Thus to completely specify the 
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FIG. 2. The tight-binding electronic band structure of bcc Mo at the LDA equilibrium lattice constant (3.12A). The energies 
of the eigenstates have been uniformly shifted to set the Fermi level to zero. 



Slater-Koster tight-binding scheme in this model requires 
the determination of at most seventy-three parameters. 
In practice, we only used the dj^ found in (|^) in the 
parametrization of Ti, Zr, and Hf, so for most elements 
we had only seventy independent parameters. The pa- 
rameters are determined by requiring that the eigenval- 
ues determined by the Slater-Koster parameteralj (0- 
reproduce the electronic band structure (including 
the eigenvalue shift (||)) produced by a set of first- 
principles density functional calculations, and that the 
corresponding energies (|^) reproduce the first-principles 
energies. For the metals discussed here we use either tha 
muffin-tin potential Augmented Plane Wave (APW)c3 
or the full|potential Linearized Augmented Plane Wave 
(LAPW)c]l'E3 method to calculate the electronic band 
structure and total energies for 4-6 volumes in each of 
the fee and bcc structures, using uniform fc-point meshes 
which include the origin and contain 89 and 55 points, 
respectively. IiV-^H cases we used the Hedin-Lundqvisi. 
parametrizatioro of the Local Density ApproximationEl 
(LDA) to DFT. Possible spin polarization is ignored. For 
a typical transition metal there are about 4000 eigenval- 
ues and energies in the database. We use an IMSL pack- 
age, based_pn a finite difference Levenberg-Marquardt 
algorithm,EJ to adjust the seventy (or seventy-three) pa- 
rameters to reproduce this database by means of a non- 
linear least squares fit, with the total energies typically 
weighted about 200 times larger than the eigenvalues in a 
single band. We place additional conditions on the mini- 
mization as follows: Since the overlap matrix defined by 
the parameters ( |l0| ) must, physically, be positive definite, 
a function which penalizes small or negative determinant 
overlap matrices is included. While this does not guaran- 
tee that our overlap matrix will remain positive definite 
during the minimization process, it does tend to keep the 



determinant from changing sign, meaning that an overlap 
matrix which is originally positive definite tends to stay 
that way. We also restrict the behavior of the functions 
(|To|). Physically, we expect both the Hamiltonian and 
Overlap parameters to decay steadily and not change sign 
for distances R near or greater than the nearest neighbor 
distance. We thus add another contribution to our min- 
imization function, which penalizes parameter sets (e-y, 
fy, g~f) which cause the either P-y{r) or P!^(r) to vanish 
over the range of distances where we are likely to require 
evaluation of the Slater-Koster matrix elements. Finally, 
we select our starting parameters guided by those found 
in previous work,t3 and ensure that the symmetry of the 
eigenstates is taken into account. 

As an illustration of the above technique, we present 
several results for molybdenum. Using the APW pro- 
gram, we calculated the total energy and band structure 
of Mo at four volumes in the fee structure and five vol- 
umes in the bcc structure. The parameters in (^, (|^), 
and ( p^ ) were then adjusted to best fit the total energies 
and the band structure in the lower six bands of Mo over 
all nine volumes. The resulting two center Slater-Koster 
parameters ( |l0[ ) are shown in Figure |l] for the Hamil- 
tonian matrix elements. We note that there is a smooth 
and monotonic decay of the TB parameters towards zero. 
To assess the quality of the fit, we calculate the tight- 
binding band structure of Mo at the LDA equilibrium 
volume (3. 12 A) of the bcc structure (Figure |^) and com^ 
pare it to the first-principles band-structure (Figure ^) .E£l 
We see that the band structures are in close agreement 
over the entire Brillouin zone and a wide spectrum of en- 
ergies, with some deterioration far above the Fermi level. 
We also calculate the energy /volume behavior for a large 
number of phases, and compare the fee and bcc ener- 
gies to those obtained from APW calculations, as shown 
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FIG. 3. The self-consistent first-principles LDA electropip band structure of bcc Mo at the LDA equilibrium lattice constant 
(3.I2A), as determined by the muffin-tin APW programE2l The energies of the eigenstates have been uniformly shifted to set 
the Fermi level to zero. 



in Figure ]^ We see that there is excellent agreement 
between the tight-binding and first-principles curves for 
the fee and bcc structures, and that aU other structures, 
including hep, are above the bcc ground state. 
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FIG. 4. The energy/volume behavior of molybdenum in 
several phases, using the tight-binding method. The energy 
of the diamond structure is not shown, as its minimum energy 
is 0.116 Ry above the origin. The small symbols (O) represent 
the APW energies of the bcc and fee phases. The construction 
of the LI2 and DO3 lattices is discussed in Table ^ The origin 
of the energy was chosen arbitrarily. 

The above procedure works well for aU of the met- 
als described here except vanadium, where the elastic 
constant C44 is predicted toJae 172 GPa, compared to 
30 GPa found experimentaUyc^ and 43 GPa when found 
using the first-principles LAPW method. This discrep- 
ancy is probably due to the high electpapic density of 
states in vanadium near the Fermi level.ta To better de- 
scribe the electronic structure of vanadium we have made 



two changes to the above formalism. First, we split the 
d on-site jenergy into contributions from the t2g and 
eg levels.Ej Second, we increase our parameter space by 
modifying the Slater-Koster parameter's functional form 
( p^ ) to include a quadratic term, 

P^{R) = (e^ + f^R + J^R^) exp[-52i?]F,(i?) . (11) 

This procedure adds an additional twenty-four param- 
eters to our basis set, for a total of ninety-seven. We 
reduce this by setting da = Q in (|^) for a = s, p, t2g^ 
and eg. When we fit the remaining ninety-three param- 
eters, as above, we reduce the predicted elastic constant 
C44 to 92 GPa without significantly changing our other 
predictions. While this is not in perfect agreement with 
experiment, it is a substantial improvement, so we will 
use these parameters for vanadium in preference to the 
seventy-three parameter set which gives a worse elastic 
constant. 

We have used this method to determine tight-binding 
parameters for most of the alkaline earth, transition and 
noble metals. In the case of the anti-ferromagnetic (Cr 
and Mn) and magnetic (Fe, Co, and Ni) elements we re- 
stricted our calculations to the non-magnetic case. This 
causes some difficulty when comparing our results to ex- 
periment, as shown below. The parameters for La are not 
presented because its / bands are located in the middle of 
the d bandsllj and therefore it is impossible to fit within 
our formalism. The Appendix lists various methods of 
obtaining the parameters used in this paper. 
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TABLE I. Equilibrium lattice constants and bulk moduli for the experimentally observed ground state structures of the 
elements, comparing the results of our tieht-binding parametrization (TB), first-principles full potential LAPW results (LAPW), 
where available, and-experiment (Exp. ).E3i23 Where not referenced, the LAPW results were done for this paper, using previously 
developed methodsB 



Element 


Structure 




a (A) 






e (A) 






Bo (GPa) 








TB 


LAPW 


Exp. 


TB 


LAPW 


Exp. 


TB 


LAPW 


Exp. 


Ca 


fee 


5.30 


5.34=^ 


5.58 


a 


a 


a 


21 


19^ 


15 


Sc 


hep 


3.26 




3.31 


4.91 




5.27 


63 




44 


Ti 


hep 


2.97 


2.87 


2.95 


4.80 


4.55 


4.68 


122 


120 


105 


V 


bee 


2.94 


2.93 


3.03 


a 


a 


a 


211 


196 


162 


Cr 


bee 


2.80'' 




2.88 


a 




a 


283 




190 


Mn 


aMn 


8.41*^ 




8.91 


a 




a 


320 




60 


Fe 


bee 


2.71'' 




2.87 


a 




a 


281 




168 


Co 


hep 


2.40'' 




2.51 


3.90'' 




4.07 


384 




191 


Ni 


fee 


3.43'' 


3.42= 


3.52 


a 


a 


a 


268 


261 


186 


Cu 


fee 


3.52 


3.52 


3.61 


a 


a 


a 


189 


190 


137 


Sr 


fee 


5.73 




6.08 


a 


a 


a 


15 




12 


Y 


hep 


3.59 




3.65 


5.35 




5.73 


46 




37 


Zr 


hep 


2.99 




3.23 


5.57 




5.15 


108 




83 


Nb 


bee 


3.25 


3.25 


3.30 


a 


a 


a 


187 


193 


170 


Mo 


bee 


3.12 


3.12= 


3.15 


a 


a 


a 


283 


291= 


272 


Tc 


hep 


2.72 




2.74 


4.34 




4.40 


304 




297 


Ru 


hep 


2.68 




2.71 


4.26 




4.28 


360 




321 


Rh 


fee 


3.77 


3.76 


3.80 


a 


a 


a 


306 


309 


270 


Pd 


fee 


3.85 


3.85 


3.89 


a 


a 


a 


212 


220 


181 


Ag 


fee 


4.01 


4.01 


4.09 


a 


a 


a 


142 


142 


101 


DBu 


bee 


4.0Z 




c no 


a 


a 


a 


1 n 
10 




1 n 
iU 


Hf 


hep 


3.07 




3.19 


5.08 




5.05 


111 




109 


Ta 


bee 


3.30 


3.24 


3.30 


a 


a 


a 


185 


224 


200 


W 


bee 


3.14 


3.14 


3.16 


a 


a 


a 


319 


333 


323 


Re 


hep 


2.78 




2.76 


4.39 




4.46 


371 




372 


Os 


hep 


2.75 




2.74 


4.31 




4.32 


441 




418 


Ir 


fee 


3.86 


3.82= 


3.84 


a 


a 


a 


389 


401= 


355 


Pt 


fee 


3.90 


3.90 


3.92 


a 


a 


a 


318 


305 


278 


An 


fee 


4.06 


4.06 


4.08 


a 


a 


a 


187 


205 


173 



^First-principles LAPW calculations 

''Non-magnetic calculation. The experimental structure exhibits some form of magnetism. 
=First-principles LAPW calculations.tl 



III. GROUND STATE BEHAVIOR AND PHASE 
STABILITY 

We first test our parameter sets by looking at the equa- 
tion of state for the observed ground state of each ele- 
ment. For the fee and bee non-magnetic cubic crystals, 
our parameters should reproduce the first-principles re- 
sults for these phases because of the fitting procedure. 
For the hep elements and manganese, these calculations 
serve as a tight-binding prediction of the equation of 
state. The case oLmanganese has been discussed in more 
detail elsewhere.E^ 

Table || shows the equilibrium lattice constants and 
bulk moduli for all of the metals studied here, Gopppar- 
ing them to fkst^principles LAPW calculationg3'El and 
to experimentEDS For the non-magnetic cubic crystals, 
the equilibrium lattice constant is within 1% of the first- 
principles LDA value for all of the cubic elements except 



tantalum and iridium, and in these cases the discrepancy 
is less than 2%. Similarly, the equilibrium bulk moduli 
are in good agreement with the first-principles results, 
within 10% for all non-magnetic cubic elements except 
tantalum, where the error is 17%. 

The hexagonal close packed phases were not fit to first- 
principles results. As a result, the equations of state pre- 
dicted by the tight-binding model are not as accurate as 
for the cubic lattices. The largest error is for zirconium, 
where a is 7% smaller than experiment and c is 8% larger. 
Yttrium and hafnium also show large discrepancies be- 
tween the tight-binding model and experiment, and the 
c for titanium is 2.5% larger than experiment. The lat- 
tice constants for the other elements are within 2% of 
experiment, consistent with the errors we would find in 
first-principles calculations. We conclude that the tight- 
binding method is nearly as good as first-principles LDA 
methods in determining the structural parameters of hep 
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TABLE II. Relative energies per atom of several structures for each of the metals examined by the tight-binding model 
discussed in the text. The energy of the experimental ground state structure is arbitrarily set to zero. All energies are 
calculated at the equilibrium volume found by the tight-binding fit, and are expressed in mRy. Below the common name of 
each phase is its Strukturbencht designation. 



Element Structure 





fee 


bcc 


hep 


diamond 


sc 


aMn 


/3Mn 


f3W 


AuCua 


AlFea 




Al 


A2 


A3 


A4 


A;, 


A12 


A13 


A15 


LI2" 


DOa"^ 


Ca 


0.0 


2.8 


0.5 


83.5 


13.3 


2.1 


3.1 


0.7 


19.8 


24.4 


Sc 


4.7 


8.5 


0.0 


101.1 


34.3 


10.2 


11.0 


19.4 


39.7 


33.9 


Ti 


5.2 


9.1 


0.0 


157.7 


64.9 


16.7 


15.4 


22.4 


56.1 


63.8 


V 


19.6 


0.0 


19.6 


180.8 


76.0 


11.2 


12.7 


12.2 


62.1 


53.4 


Cr 


28.6 


0.0 


30.6 


49.1 


118.6 


20.3 


19.7 


16.2 


94.3 


87.2 


Mn 


7.5 


14.3 


2.9 


44.2 


90.3 


0.0 


1.5 


7.0 


87.6 


66.2 


Fe 


-27.9 


0.0 


-35.0 


79.8 


54.5 


-23.1 


-25.0 


-8.1 


4.5 


20.4 


Co 


-2.9 


21.0 


0.0 


135.0 


8.2 






19.2 


72.7 


75.2 


Ni 


0.0 


8.0 


2.8 


142.3 


75.9 


7.7 


4.6 


15.0 


48.7 


57.1 


Cu 


0.0 


3.5 


1.2 


70.8 


24.7 


6.2 


5.4 


10.1 


14.7 


22.0 


Sr 


0.0 


2.1 


1.0 


73.0 


22.0 


3.2 


3.2 


1.4 


22.2 


25.4 


Y 


2.4 


8.8 


0.0 


97.9 


432.2 


8.6 


10.4 


17.3 


31.5 


27.5 


Zr 


1.1 


2.9 


0.0 


112.7 


123.3 


8.6 


8.0 


8.0 


18.7 


26.2 


Nb 


26.1 


0.0 


22.7 


181.9 


68.6 


17.1 


21.1 


15.9 


69.5 


59.6 


Mo 


30.0 


0.0 


31.0 


147.3 


68.7 


17.5 


19.9 


11.3 


67.7 


53.5 


Tc 


6.1 


23.3 


0.0 


72.9 


56.3 


0.2 


4.5 


13.3 


64.6 


35.4 


Ru 


1.0 


44.6 


0.0 


134.5 


106.6 


16.3 


14.6 


34.7 


103.0 


82.7 


Rh 


0.0 


31.7 


5.1 


159.8 


96.0 


16.3 


7.4 


29.5 


81.2 


88.3 


Pd 


0.0 


9.9 


2.5 


148.6 


84.9 


9.2 


6.1 


12.8 


50.1 


72.8 


Ag 


0.0 


3.3 


0.6 


66.1 


24.1 


6.0 


5.3 


8.6 


15.7 


25.2 


Ba 


0.9 


0.0 


0.1 


85.3 


23.6 


2.8 


2.9 


2.2 


27.0 


30.7 


Hf 


1.2 


7.1 


0.0 


380.2 


113 


22.9 


20.6 


20.8 


110.2 


143.3 


Ta 


24.7 


0.0 


25.3 


188.6 


63.2 


8.9 


14.0 


9.8 


65.6 


54.2 


W 


36.0 


0.0 


10.1 


159.3 


113.7 


18.2 


25.7 


12.4 


134.2 


94.2 


Re 


13.4 


27.7 


0.0 


28.6 


55.6 


0.2 


1.8 


19.1 


61.0 


22.1 


Os 


8.1 


64.3 


0.0 


11.2 


58.6 


17.8 


17.4 


42.4 


66.2 


33.6 


Ir 


0.0 


50.1 


8.5 


83.9 


83.5 


23.8 


19.9 


41.8 


64.0 


75.6 


Pt 


0.0 


9.7 


4.4 


168.7 


78.3 


14.4 


12.8 


25.9 


48.5 


78.5 


An 


0.0 


1.4 


0.6 


96.1 


24.1 


8.7 


8.4 


15.2 


17.5 


32.3 



^Calculations using the LI2 (AuCua) lattice, placing the metal atoms on the copper sites and leaving vacancies on the gold 
sites. 

''Calculations using the DOa (AlFea) lattice, placing the metal atoms on the iron sites and leaving vacancies on the aluminum 
sites. 



metals. 

To be successful, the tight-binding method must also 
show that the experimental ground state structure is, 
in fact, the ground state predicted by the method. We 
have calculated energy versus volume curves for several 
common crystal structures. As an example, Figure ^ 
shows the resulting energy/ volume curves for molybde- 
num. The results for all of the elements are summa- 
rized in Table ||, which shows the equilibrium energy of 
each of these phases expressed as the difference in en- 
ergy between that phase and the equilibrium energy of 
the experimental ground state. The c/a ratio of the hep 
phase has been chosen to minimize the tight-binding en- 
ergy at each volume to obtain the hep energies. The 
tight-binding method correctly predicts the ground state 
structure for all metals except the ferromagnetic metals 
iron and cobalt. This is understandable, since we have 



not included spin-polarization. Especially satisfying is 
the fact the method finds the correct ground states of 
the hep metals and manganese,E3 since these phases are 
not included in the fit of the tight-binding parameters. 

We can also assess the stability of the tight-binding pa- 
rameters by studying crystal structures of lower symme- 
try. Thus Figure ^ shows the energy of molybdenum un- 
der a volupie conserving tetragonal strain, the so-called 
Bain pathp3 at the experimentally observed equilibrium 
volume of the bcc phase. The energy is properly a min- 
imum for the bcc lattice, where c/a — 1, and attains 
a local maximum at the fee structure (c/a = V2), in- 
dicating that Cii — C12 < 0, so that the fee structure 
is unstable. Similarly, Figure ^ shows the energy for a 
monatomic lattice undergoing a trigonal strain with the 
primitive vectors 
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FIG. 5. The tight-binding calculation of the energy of 
molybdenum, at the experimental equilibrium volume, un- 
der a tetragonal strain (the Bain pathc3), as a function of 
c/a. The vertical lines denote the positions of the bcc and fee 
lattices. 



ai \ I a 13 13 
as / \pf3a 




(12) 



where a and (3 are chosen so that the volume is the ex- 
perimentally observed volume of the bcc lattice, and the 
angle between the primitive vectors is given by 9, where 



cos 9 — 



(2a + I3)p 
2/?2 



(13) 



The vectors (g^ represent an fee lattice when 9 = 60 deg, 
a simple cubic lattice when 9 = 90 deg, and a bcc lattice 
when 9 = 109.47 deg. The tight-binding calculations in 
Figure |6| show that the bcc phase is, indeed, the minimum 
energy structure. Also, the simple cubic phase has a local 
maximum in the energy, indicating that it is elastically 
unstable with C44 < 0. 



IV. ELASTIC CONSTANTS 

Elastic constants are calculated by imposing an exter- 
nal strain on the crystal, relaxing any internal parame- 
ters to obtain the energy as a function of the strain, and 
numerically solving for the elastic constaats as the cur- 
vature of the energy versus strain curve.cTEl We have cal- 
culated the elastic constants for the ground state phases 
of all of the elements in this study except manganese. 
This procedure serves two purposes: it tests the mechan- 
ical stability of the tight-binding ground state; and it as- 
sesses the ability of the method to determine properties 
not included in the fit of the tight-binding parameters. 

The elastic constants of the cubic materials (except 
manganese) are shown in Table 111. This table also 
shows the elastic constants obtained frnrn. fii'|St-principles 
LAPW computationsQ and experiment .E^Ej'Ell All calcu- 
lations were performed at the experimental volume. The 




80 90 
e (deg) 

FIG. 6. The tight-binding calculation of the energy of 
molybdenum, at the experimental equilibrium volume, un- 
der a trigonal strain of the lattice given by (^2|) , as a function 
of the angle between the primitive vectors, 9. The positions 
of the bcc, simple cubic, and fee lattices are denoted by the 
vertical lines. 



accuracy of the elastic constants Cn and C12 is remark- 
ably good. For the twelve elements for which we have 
first-principles results, none of the tight-binding elastic 
constants differs from the first-principles results by more 
than 15%. For comparison, the first-principles results 
also differ from the experimental results by less than 
15%, a measure of the accuracy of the LDA in determin- 
ing these elastic constants. Thus the tight-binding elas- 
tic constants are not significantly worse than those com- 
puted from first-principles. Considering only the non- 
magnetic materials, the root-mean-square (RMS) devia- 
tion of the tight- binding Cn compared to experiment is 
38 GPa, while the RMS deviation of C12 is 15 GPa. 

The errors in C44 are somewhat larger, the RMS devia- 
tion from experiment being 49 GPa. In this case, if we ex- 
clude calcium and strontium, discussed below, the maxi- 
mum deviation of first-principles results from experiment 
is 22%, in molybdenum. The maximum deviation of the 
tight-binding clastic constant from the first-principles re- 
sult is in vanadium, where the error is 114%, even when 
we use the quadratic Slater-Koster parametrization ( |Tl| ) 
and split the on-site d terms. As noted above, this is 
probably due to the high value Q^-ithe density of states of 
vanadium near the Fermi levelJlj Excluding vanadium, 
the maximum deviation is 36%, in niobium. For compar- 
ison, the first-principles results differ from experiment 
by at most 22%, in molybdenum. Hence, except for the 
anomaly of vanadium, the tight binding results are only 
slightly worse than the first-principles results in the de- 
termination of elastic constants. 

Calculation of C44 in the alkaline earth metals is dif- 
ficult because the lattice is extremely soft. The error in 
calcium is 71%, compared to first-principles results, and 
50% compared to experiment. For strontium the situa- 
tion is even worse, since we predict the fee lattice to be 
unstable to any strain which is related to C44. In this 
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TABLE III. Elastic constants for cubic elements. All ele- 
ments for which we have constructed a set of tight-binding 
parameters and which have a cubic ground state (except 
manganese) are presented here. Comparison is made be- 
tween the results of our tight-binding parametrization (TB), 
first-principles full potential LABW results (LAPW), where 
available, and experiment (Exp.)£j Where unreferenced, the 
LAPW results ware done for this paper, using previously de- 
veloped methods.u Calculations were performed at the exper- 
imental volume. 





TB 


oil 
LAPW 


Exp 


TB 


LAPW 


Exp 


TB 


LAPW 


Exp 




10 


1 / 


ID 


lU 


lU 


IZ 




14 


o 


V 




OOP; 




lUD 


111 
111 


1 1 o 

iiy 


OO 


Qn 
oU 


A Q 
4o 


or 








oo 




DD 


zou 




lUU 


Fe 


-4 






227 






180 






Ni 


256 






142 






86 






Cu 


161 


156 


168 


108 


106 


121 


55 


80 


75 


Sr 


8 




15 


3 




6 


-3 




10 


Nb 


204 


230 


246 


137 


122 


139 


34 


25 


29 


Mo 


453 


468^ 


450== 


147 


149'' 


173" 


120 


98=" 


125^^ 


Rh 


491 


433 


413 


171 


185 


194 


260 


206 


184 


Pd 


233 


218 


227 


163 


172 


176 


63 


74 


72 


Ag 


133 


122 


124 


86 


90 


93 


42 


52 


46 


Ba 


9 






7 






13 






Ta 


275 


256 


261 


140 


154 


157 


78 


67 


82 


W 


529 


527 


523 


170 


194 


203 


198 


147 


160 


Ir 


694 


62r 


590 


260 


256^ 


249 


348 


260^ 


262 


Pt 


380 


381 


347 


257 


189 


251 


71 


83 


76 


Au 


184 


200 


189 


154 


173 


159 


43 


33 


42 



^First-princif 
''Experiment 
"^ExperimentE 



Ks LAPW calculations.! 



regard, however, it is important to note that even many 
of the best LDA elastic constant calculations in Table [II 
have an absolute error of more than 10 GPa. Since the 
experimental values of C44 for both calcium and stron- 
tium are less than 10 GPa, it is understandable that we 
cannot reproduce the experimental elastic constants for 
the alkaline earth metals with any degree of accuracy. 

We also note that this method correctly reproduces the 
sign of the elastic constant C12 — C44, even in rhodium 
and iridium, where it is negative, though the predicted 
magnitude \Ci2 — C44I is much larger than observed ex- 
perimentally and in first-principles calculations for these 
two materials. The negative sign cannot be obtained 
from the standard Embedded Atom Method (EAM),l3 
though it can be faujid in the Modified Embedded Atom 
Method (MEAM).EI 

We have also calculated elastic constants for all of the 
non-magnetic transition metals which take the hep struc- 
ture. Hexagonal close packed crystals have two atoms 
per unit cell, and the internal parameter describing the 
atomic positions is free to move when we strain the crys- 
tal to calculate elastic constants. Using the tight-binding 
approach, it is relatively simple to determine the value 



TABLE IV. Elastic constants for hexagonal close-packed 
elements. All elements for which we have constructed a set 
of tight-binding parameters and which have a non-magnetic 
hep ground state are presented here. Comparison is made be- 
tween the results of oipi tight-binding parametrization (TB) 
and experiment (Exp) SB The tight-binding results include in- 
ternal relaxation. Calculations were performed at the exper- 
imental volume, but at the c/a ratio which minimized the 
energy for that volume. Note that in a hexagonal crystal, 
Cee = {Cii — Ci2)/2. 





Cii 


C12 


C13 


C33 


C44 




TB 


Exp 


TB 


Exp 


TB 


Exp 


TB 


Exp 


TB 


Exp 


Sc 


73 


99 


30 


40 


31 


29 


78 


107 


25 


28 


Ti 


171 


160 


58 


90 


46 


60 


203 


181 


64 


47 


Y 


45 


78 


15 


29 


13 


20 


48 


77 


17 


24 


Zr 


102 


144 


67 


74 


65 


67 


99 


166 


4 


33 


Tc 


477 




196 




172 




505 




127 




Ru 


642 


563 


170 


188 


129 


168 


706 


624 


214 


181 


Re 


559 


616 


283 


273 


250 


206 


621 


683 


136 


161 


Os 


754 




272 




274 




831 




232 





of the internal parameter which minimizes the total en- 
ergy for a given external strain, and then find the related 
elastic constant by standard techniques. This is more 
difficult for first principles calculations, so we have not 
calculated LAPW elastic constants for hexagonal materi- 
als. The tisl^t-binding elastic constants are compared to 
experimenttj in Table For simplicity in interpreting 
the results, we calculated the elastic constants using the 
tight-binding c/a ratio at the experimental volume. The 
hep tight-binding elastic constants show larger relative 
deviations from experiment than found in the cubic crys- 
tals. The largest deviation, 88%, occurs in the calculation 
of C44 in zirconium, where we predict a value of 4 GPa, 
compared to the 33 GPa found by experiment. The RMS 
deviations of the tight-binding elastic constants from ex- 
periment are 46 GPa for Cn, 17 GPa for C12, 25 GPa 
for Ci3, 54 GPa for C33, and 22 GPa for C44. The RMS 
deviation over all of the hep elastic constants is 36 GPa. 
When compared to the RMS deviation for cubic mate- 
rials (37 GPa), we see that on average the hep elastic 
constants are as accurate as those calculated in the cubic 
case. 



V. VACANCIES 

Vacancy formation energies are mo^^easily determined 
by the supercell total-energy method. a One atom in the 
supercell is removed, and neighboring atoms are allowed 
to relax around this vacancy while preserving the sym- 
metry of the lattice. The great advantage of the tight- 
binding method over first-principles calculations is that 
we can do the calculation in a very large supercell, elim- 
inating the possibility of vacancy-vacancy interactions. 
We find that a supercell containing 108 atoms is sufR- 
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TABLE V. Tight-binding vacancy formation energies com- 
pared to first-principles calculations and experiment. Ener- 
gies were computed using a 108 atom supercell. Calculations 
with the atoms at the primitive lattice sites in the crystal 
(Fixed) and allowing relaxation around the vacancy (Relaxed) 
are shown. First-principles calculations of the vacancy forma- 
tion energy are given in the column labeled "LDA". The ex- 
perimental column shows a range of energies if several experi- 
ments have been tabulated. Otherwise the estimated error in 
the experiment is given. 

Vacancy Formation Energy 



Element 


Tight-Binding 


LDA 


Experiment" 




Fixed 


Relaxed 






Cu 


1.29 


1.18 


1.41'" 1.29"^ 


1.28 - 1.42 


Nb 


2.84 


2.82 




2.65 ± 0.3 


Mo 


2.63 


2.46 




3.0 - 3.6 


Rh 


3.39 


3.35 


2.26"^ 


1.71"^ 


Pd 


2.46 


2.45 


1.57*= 


1.85 ± 0.25 


Ag 


1.31 


1.24 


1.20*^ 1.06'' 


1.11 - 1.31 


Ta 


3.17 


2.95 




2.9 ± 0.4 


W 


6.86 


6.43 




4.6 ± 0.8 


Ir 


2.19 


2.17 




1.97'' 


Pt 


2.79 


2.79 




1.35 ± 0.09 


Au 


1.24 


1.12 




0.89 ± 0.04 



"Experimental energies as tabulated by Schaefe: 
erwise noted. 

''First-principles calculations of Dederichs et ai. 
"First-principles calculations of Polatoglou et al. 
''Experimental energy tabulated by de Boer et al 



unless oth- 



cient to eliminate the interaction. The vacancy formation 
energy is given by 

EyaciV) = EsciN - 1, 1; V)-{N- l)Ei,uik{V/N) , 

(14) 

where Esc{M, Q; V) is the total energy of a supercell of 
volume V containing M atoms and Q vacancies, and 
EbuikiV) is the energy per atom of the bulk metal at 
a volume V per atom. We use the experimental lattice 
constant to set the volume of the system, since under 
experimental conditions the lattice constant of a metal 
containing isolated vacancies will be the lattice constant 
of the bulk metal. 

We have computed the vacancy formation energy of 
several of the cubic transition and noble metals using 
this tight binding method, including relaxation around 
the vacancy. The results are presented in "feb^ 0, 
where we also jCompare to first principles resultg23't3 and 
experiment EJ'Ej The calculated vacancy formation ener- 
gies of niobium, silver, tantalum, and iridium are in ex- 
cellent agreement with experiment. The tight-binding 
parametrizations of rhodium, tungsten, and platinum 
give the poorest vacancy formation energies, 1-2 eV 
above the experimental value. 

We note that the LI2 and -DO3 lattices described in 
Table |l| can be considered vacancy-containing supercells 



with N = A. Though the vacancy-vacancy interaction in 
these cells is quite large, fitting the TB parameters ([7[]lO|) 
to first-principles results for these lattices will improve 
the vacancy formation energy calculation. We have yet 
to implement this procedure. 



VI. SURFACES 

The tight-binding model can be used to calculate sur- 
face energies by the supercell technique. A slab of metal 
is formed by cleaving the crystal along the desired plane, 
creating two identical free surfaces. The distance between 
the two surfaces in increased, creating a set of slabs which 
repeat periodically in the direction perpendicular to the 
surfaces. Of course, we must take care to separate the 
slabs by a large region of vacuum so that the electrons 
on one slab cannot hop to a neighboring slab. The slabs 
must be thick enough so that the atoms at the center of 
the slab have the properties (e.g., local electronic density 
of states) of atoms in the bulk material, and so that the 
two surfaces on the same slab cannot interact with each 
other. For low index ((100), (110), and (111)) faces of 
the non-magnetic fee and bcc metals we find that a slab 
containing 25 atomic layers is sufficient to meet these cri- 
teria. We also converge the calculations with respect to 
the /c-point mesh. Depending on the surface and underly- 
ing bulk structure this takes between 19 and 91 /c-points 
in the two dimensional Brillouin zone. We estimate nu- 
merical error in the TB surface energies to be about 0.1 
J/m2. 

The surface energy, expressed as the energy required 
to create a unit area of new surface, is then given by the 
formula 



Esurf — -^{Eslab — N E^ulk) 



(15) 



where A is the area occupied by one unit cell on the 
surface of the slab, Esiab is the total energy of the slab, 
N is the number of atoms in the unit cell, and E^uik is 
the energy of one atom in the bulk at the lattice constant 
of the atoms in the interior of the slab. For simplicity, 
the calculations were done using the bulk equilibrium 
lattice parameters, with no relaxation or reconstruction 
at the surface. Our calculations indicate that relaxation 
energies are on the order of 1-10% of the total surface 
energy, so neglecting relaxation will not significantly alter 
our conclusions. Of course, relaxation and reconstruction 
are necessary in order to understand surface properties. 
We will address these issues in a future paper. 

We have calculated the formation energies for bulk ter- 
minated (i.e., unrelaxed and unreconstructed) low index 
((100), (110), and (111)) surfaces of the non-magnetic 
fee and bcc transition and noble Hietals. Table VI 
compares these rppilts to experiment first-nriTjiciples 
calculations ,L3oCj and atomistic models.L3c30 Note 
that except for the (100) surface of tungsten, the experi- 
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TABLE VI. Surface energies, calculated from the tight-binding theory (TB) , by first principles Local Density methods (LDA) , 
by the embedded atom method (EAM), or by the modified embedded atom method (MEAM), compared to experiment. Energies 
are given in units of J/m? . 



Element 


Orientation 


TB 


LDA 


EAM 


MEAM" 


Experiment 




(111) 


3.34 






1.81 


2.6'= 


V 


(110) 


2.13 




1.68^^ 


1.71 






(100) 


3.04 




1.83" 


2.49 






(111) 


1.73 


1.94'' 


1.17^^ 


1.41 


1.77' 


Cu 


(110) 


2.04 




1.40" 


1.64 






(100) 


1.93 




1.28° 


1.65 






(111) 


2.44 






2.02 


2.3'= 


Nb 


(110) 


1.54 


2.36*5 2.9'' 


1.8^ 


1.87 






(100) 


2.37 


2.86*5 3.1'' 


1.97" 


2.02 






(111) 


2.84 






1.86 


2.9" 


Mo 


(110) 


3.04 


3.14*5 


2.13" 


1.93 






(100) 


2.12 


3.52*5 


2.28" 


2.12 






(111) 


2.46 


2.54" 2.53*5 




2.60 


2.6" 


Rh 


(110) 


2.71 


2.88*5 




2.92 






(100) 


2.57 


2.81*5 




2.90 






(111) 


1.57 


1.64*5 


1.22" 


1.38 


2.00" 


Pd 


(110) 


1.86 


1.97*5 2.5'' 2.70' 


1.49" 


1.67 






(100) 


1.75 


1.86*5 2.3'' 1.04' 


1.37" 


1.66 






(111) 


1.14 


1.21*5 121" 


0.62" 


1.09 


1.32' 


Ag 


(110) 


1.42 


1.26*5 1.4'' 


0.77" 


1.27 






(100) 


1.29 


1.21*5 ;13h 


0.71" 


1.22 






(111) 


3.14 






2.31 


2.78" 


Ta 


(110) 


2.05 




1.80" 


2.17 






(100) 


3.00 




1.99" 


3.29 






(111) 


6.75 






2.25 


2.99' 


W 


(110) 


4.30 




2.60" 


2.23 






(100) 


6.7 


5.54j 


2.81" 


2.65 


6'^ 




(111) 


2.59 






2.84 


3.0" 


Ir 


(110) 
(100) 


3.19 

2.95 






3.06 
2.91 






(111) 


2.51 




1.44" 


1.66 


2.49" 


Pt 


(110) 


2.97 




1.75" 


2.13 






(100) 


2.83 




1.65" 


2.17 






(111) 


1.48 




0.79" 


0.89 


1.54' 


Au 


(110) 


1.85 




0.98" 


1.12 






(100) 


1.69 




0.92" 


1.08 





"Modified Embedded- Atom calculationsEH . 
'^Experimental determination of the surface eBergy of an "average" face using an approximate extrapolation to T = 0.1 
"Embedded- Atom calculations for bcc metals.tJ |— ■ 
"Full-Potential Linearized Muffin- Tin OrbitaLiFuU Potential-LMTO) calculation, using seven layer slabs.E3 
"Embedded- Atom calculations for fee metals.LJ I — . 
'Experimental determination of the surface eneifgt' of an "average" face at T = OO 
^FP-LMTO calculation, using seven layer slabs.eS |— . 
''Full-Potential LAPW calculation, using nine layer slabs.O i— i 
'Local orbital pseudo-potential calculation, using a three layer slab.Ej 

^Full-Potential Linearized Augmented PlarpiWave (LAPW) calculation for the ideal surface, using a five layer slab 
''Experimental result for the (100) surfacctS 
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mental resultg23'Ea are for an "average" surface, and|Some 
are extrapolated to T = from high temper aturesS 

The results for the fee metals are particularly gratify- 
ing. In all cases we find < i?(ioo) < -E'(iio)- Thus 
close-packed surfaces are the most stable for the fee met- 
als. The TB surface energies are uniformly larger, and 
closer to experiment, thaij-those obtained by the embed- 
ded atom method (EAM),c3 which is known to underesti- 
mate surface energies in fee metals. Our surface energies 
are generally closer to experiment than those obtained 
by the Modified Embedded Atom Method (MEAM),y 
and in good agreement with first-priuciples calculations, 
of the surface energy of copper (111)J23 and rhodium.c3 
Our calculations are also in the range of the surface en- 
ergies predicted by sesieijal first-principles calculations for 
palladium and silver .cd^cS 

Results for bcc metals are also good. For all of the 
unrelaxed surfaces we find -E(iio) < £-(100) < -E^(iii)- The 
energies for niobium and molybdenum are somawjhat low 
compared to both first-principles calculationsEdcj and 
experiment, but the only serious discrepancy is the aver- 
age surface energy of tungsten, which is smaller than our 
values by a factor of one-third. At the moment this is not 
particularly worrisome, because experimental estimates 
of the surfacej-energy of tungsten range from 1.68 J/m^ 
to 4.50 J/m^O The latter number is, in fact, in range 
of our calculation for the (110) surface. In addition, the 
results for the energy of the (100) surface of tungstea 
are in good agreement with both first-principles resultsufl 
and experiment .Lj. In most cases the tight-binding sur- 
face energies are greater than either of the corresponding 
EAM or MEAM values. Where this is not the case the 
energy difference is slight. 



Compared to experiment, the accuracy of the vacancy 
formation energy predictions range from good to poor, 
with niobium, silver, tantalum, and iridium all in excel- 
lent agreement with experiment, and rhodium, tungsten, 
and platinum much larger than experiment. 

Surface energies for both fee and bcc crystals are in 
good agreement with the available first-principles and ex- 
perimental data. Energies for fee crystals are uniformly 
larger than the corresponding EAM energies, in agree- 
ment with experiment. 

Improvements to the tight-binding method are un- 
der investigation. In particular, we are examining other 
forms for the parametrization of the on-site terms (|^) 
and the Slater-Koster terms (p^. We are also investigat- 
ing the effect of two-center Slater-Koster integrals-pn the 
on-site terms, as suggested by Mercer and Chou.ta Since 
these terms properly account for the crystal-field split- 
ting of the on-site terms, they may be useful in obtaining 
even more accurate results. Finally, we have begun the 
process of. extending the tight-binding method to binary 
systems. Ej 
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APPENDIX: THE TIGHT-BINDING 
PARAMETERS 



VII. SUMMARY 

We have extended the tight-binding methodS to in- 
clude most of the alkaline earth, transition, and noble 
metals. Tight-binding parameters are determined by do- 
ing a least squares fit to simultaneously reproduce the 
total energy and electronic structure determined by first- 
principles at several volumes of the fee and bcc lattices. 
In all cases except the ferromagnetic metals the method 
correctly predicts the ground state structure, even when 
the ground state is hep or aMn. 

The method was tested against first-principles results 
and experiment by the calculation of elastic constants, 
vacancy formation energies, and surface energies. Elastic 
constants for non-magnetic cubic materials are, in gen- 
eral, in good agreement with both LAPW calculations 
and experiment, though we have some problems with C44 
in vanadium. Since the tight-binding method does not 
require a specific sign for C12 — C44 we get reasonable 
elastic constants even for rhodium and iridium. Elas- 
tic constants for the hep phase are nearly as accurate as 
those calculated in the cubic phase. 



The tight-binding parameters used in this paper are 
available from the AIP Physics Auxiliary Publica tion 



Service, tittp : //www . aip . org/ epaps/ epaps .html 



as 



no ted in a footnote to this paper, on the W orld Wide Web 
at http : / / cst-www .nrl .navy .mil/bind, or via email at 
mehlSdave . nrl . navy . mil. 
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TABLE VII. Tight-binding parameters for Calcium, Scandium, Titanium, Vanadium, and Chromium. The equation 

numbers refer to the parametrization forms discussed in the text. Note that the t2g and eg on-site parameters are identical for 
all elements except vanadium. 



Parameter 


Calcium 


Scandium 


Titanium 


Vanadium 


Chromium 


On-Site Parameters Eqs. (?|) and (3) 


A 


1.12775E+00 


1.16238E+00 


1.30248E+00 


1.578632+00 


1.59665E+00 


as 


2.43912E-01 


3.91870E-01 


-1.53116E-02 


1.24744E-01 


3.57154E-01 


bs 


-1.73150E-h01 


-1.29635E+01 


1.82222E+01 


-3.46195E+01 


-8.83822E+00 


Cs 


5.96806E-F02 


4.06795E+02 


2.12586E+03 


6.62479E+04 


3.91628E+04 


ds 


O.OOOOOE-FOO 


O.OOOOOE+00 


3.35194E+03 


O.OOOOOE+OO 


O.OOOOOE+OO 


Up 


1.67761E-01 


1.08349E-01 


3.73246E-01 


5.09467E-01 


6.65273E-01 


bp 


1.87785E-h01 


1.49373E+01 


2.87173E+01 


1.22445E+02 


9.54968E+01 


Cp 


-2.11382E-f02 


-6.25612E+00 


-6.66496E+01 


-5.81167E+03 


4.05266E+02 


dp 


O.OOOOOE-HOO 


O.OOOOOE+OO 


-1.89929E+03 


O.OOOOOE+OO 


O.OOOOOE+00 


at2g 


2.31793E-01 


1.92461E-01 


1.44001E-01 


1.35087E-01 


9.80248E-02 


bt2a 


2.56001E-h00 


-4.09038E-02 


5.62900E+00 


2.82780E+01 


3.16359E+01 




4.43144E-f00 


4.49403E+01 


-1.11934E+02 


1.96581E+02 


-1.70206E+03 


dt2g 


O.OOOOOE-hOO 


O.OOOOOE+OO 


4.39451E+03 


O.OOOOOE+OO 


O.OOOOOE+OO 


Cleg 


2.31793E-01 


1.92461E-01 


1.44001E-01 


1.52188E-01 


9.80248E-02 


beg 


2.56001E-F00 


-4.09038E-02 


5.62900E+00 


1.42941E+01 


3.16359E+01 


Ceg 


4.43144E-H00 


4.49403E+01 


-1.11934E+02 


2.92193E+03 


-1.70206E+03 


deg 


O.OOOOOE+00 


O.OOOOOE+OO 


4.39451E+03 


O.OOOOOE+OO 


O.OOOOOE+OO 






Hamiltonian Parameters Eq. (lOl) 


or (11) 






1.03105E-f02 


8.97124E+01 


-3.04758E+01 


3.36262E+00 


3.16386E+00 


fsscr 


-1.75892E-K01 


-2.00737E+01 


3.65365E-01 


-1.14810E+00 


-1.12865E+00 


7 

J SS(T 


O.OOOOOE-FOO 


O.OOOOOE+OO 


O.OOOOOE+OO 


-7.17525E-03 


O.OOOOOE+00 


QsSCT 


9.40409E-01 


1.03023E+00 


1.02322E+00 


8.06897E-01 


8.40762E-01 


^SpCT 


-6.67410E-h00 


-5.39161E+00 


2.52625E+00 


9.10588E-01 


1.42516E+00 


fspa 


1.23384E+00 


1.28676E+00 


-6.91609E-02 


-5.48139E-02 


-5.72265E-02 


f spa 


O.OOOOOE-HOO 


O.OOOOOE+OO 


O.OOOOOE+00 


2.14256E-02 


O.OOOOOE+00 


Qspa 


7.43630E-01 


7.90720E-01 


7.98877E-01 


7.28117E-01 


6.78224E-01 


Gppa 


-2.74551E-01 


-1.23867E-01 


1.15348E+00 


-5.58852E-01 


1.68736E+00 


fppo- 


1.31541E-01 


1.08003E-01 


-6.14491E-02 


1.50760E-01 


1.15999E-02 


f ppa 


O.OOOOOE+00 


O.OOOOOE+OO 


O.OOOOOE+OO 


4.95827E-01 


O.OOOOOE+00 


Qppa- 


6.19331E-01 


6.19587E-01 


5.59035E-01 


7.59313E-01 


7.60586E-01 


GppTT 


-4.29276E-H02 


-2.38965E+01 


-3.89463E+00 


2.68584E+02 


-4.81928E+01 


fppiT 


4.41085E-H01 


-8.63951E+00 


1.66898E+00 


-5.61785E+01 


1.44596E+01 


7 

■/ ppn 


OOOOOE+00 


O.OOOOOE+OO 


O.OOOOOE+OO 


1 24700E+00 


OOOOOE+00 


Qppn 


1.16964E-F00 


1.14561E+00 


1.04814E+00 


1.26040E+00 


1.29469E+00 




9.68472E-01 


5.15292E-01 


-1.07381E+00 


1.41944E-02 


7.71274E-01 


fsda 


-2.37815E-01 


-1.44198E-01 


6.14641E-02 


-5.56681E-02 


-1.62549E-01 


f sda 


O.OOOOOE+00 


O.OOOOOE+OO 


O.OOOOOE+OO 


-7.17105E-04 


O.OOOOOE+00 


gsdcr 


7.16040E-01 


7.61250E-01 


7.50489E-01 


6.69904E-01 


7.41397E-01 


dpda 


1.26136E+00 


4.44491E-01 


-3.22172E-01 


4.32047E+00 


4.28533E+00 


fpdfT 


-2.09090E-01 


-1.06993E-01 


-2.21047E-03 


-1.09109E+00 


-1.06533E+00 


f pdfj 


Q.OOOOOE+00 


O.OOOOOE+OO 


O.OOOOOE+00 


-3.56750E-03 


O.OOOOOE+00 


Qpda 


6.70733E-01 


6.36292E-01 


6.24299E-01 


8.19499E-01 


8.40053E-01 


C^pdiT 


5.57017E-01 


5.24971E-01 


3.19996E-01 


-2.39826E+00 


-2.29533E+00 


fpdn 


-4.67466E-02 


-2.50914E-02 


-2.06480E-02 


7.87639E-01 


7.92622E-01 


f pd-iT 


O.OOOOOE+00 


O.OOOOOE+OO 


O.OOOOOE+OO 


-3.28881E-03 


O.OOOOOE+00 


Qpd-R 


6.00789E-01 


6.55945E-01 


6.50894E-01 


8.95260E-01 


8.52608E-01 
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4.76025E+00 


9.13632E-01 


-8.83254E-01 


-1.38932E+00 


-1.44185E+00 


fdda 


-9.11415E-01 


-2.84308E-01 


-1.06276E-02 


4.40043E-02 


8.41726E-02 


f ddtj 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+00 


-9.20172E-03 


O.OOOOOE+00 


Qdda 


8.17314E-01 


7.58386E-01 


7.48868E-01 


8.12780E-01 


7.92347E-01 




-8.83124E-01 


8.09225E-01 


2.14981E+00 


2.44958E+00 


3.37805E+00 


fddn 


5.71572E-01 


3.86043E-01 


4.84974E-02 


1.90258E-01 


3.41379E-02 


f ddiT 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+OO 


1.08500E+00 


O.OOOOOE+OO 


QddTT 


8.35313E-01 


8.67158E-01 


8.78521E-01 


9.45587E-01 


9.58150E-01 


GddS 


1.05244E+02 


1.92969E+02 


1.26999E+01 


1.52551E+02 


-1.14015E+01 


fddS 


-1.89498E+01 


-3.97738E+01 


-1.04051E+01 


-1.85124E+02 


5.73356E-01 


fddS 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+OO 


-1.58858E+00 


O.OOOOOE+OO 


QddS 


1.07735E+00 


1.18457E+00 


1.27903E+00 


1.52346E+00 


1.18366E+00 






Overlap Parameters Eq. (LO) 


or (11) 




Gssa 


7.92707E+00 


4.81308E+00 


1.13429E+00 


-6.01744E-01 


-1.90233E+00 


fssa 


-7.22811E-01 


-5.79825E-01 


-1.19576E-02 


9.74987E-01 


1.28953E+00 


7 

■1 ssa 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+OO 


2.92855E-03 


O.OOOOOE+00 


Qssa 


6.64862E-01 


7.13371E-01 


6.96534E-01 


7.63881E-01 


8.44029E-01 




5.42081E+04 


2.63307E+04 


-6.19606E+01 


9.22403E-01 


6.30159E+00 


fspcr 


-8.78337E+03 


-5.56247E+03 


5.18187E-01 


-1.85853E+00 


-1.59725E+00 


7 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+OO 


-2.50567E-01 


O.OOOOOE+00 


QaptT 


1.22836E+00 


1.35254E+00 


1.13620E+00 


8.16578E-01 


8.56016E-01 


(^ppa 


2.93593E+01 


6.97663E+00 


-3.75786E+01 


1.27358E+01 


8.43048E+00 


fppa 


-6.62151E+00 


-1.72009E+00 


1.08078E+01 


-3.00275E+00 


-2.13782E+00 


f 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+OO 


-2.76667E-02 


O.OOOOOE+00 


QpptT 


8.46628E-01 


7.93212E-01 


1.04440E+00 


7.60395E-01 


7.01353E-01 


^ppTT 


-2.47333E+02 


-1.43227E+01 


4.34031E+00 


-2.56081E+00 


-2.19353E+00 


fpp-IT 


4.47028E+01 


3.41589E+00 


-3.75328E-03 


1.61177E+00 


1.07267E+00 


f ppiT 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+OO 


-1.07454E-01 


O.OOOOOE+00 


Qpp-n 


9.65850E-01 


8.08203E-01 


8.60463E-01 


7.30749E-01 


7.56032E-01 




-1.38111E+02 


3.70841E+00 


8.62354E-01 


-1.97293E+02 


-6.35844E+02 


fsda 


2.56352E+01 


-2.31465E-01 


1.49029E-01 


4.13772E+01 


1.54736E+02 


f sda 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+OO 


7.17648E-01 


O.OOOOOE+00 


QsdtT 


9.15974E-01 


fi 9,Sfil,'iE— 01 


7 f)3941 K— ni 


1.09374E+00 


1.24363E+00 


^pda 


-7.63060E+01 


-3.42076E+02 


-1.98379E+00 


-1.06522E+01 


-1.66486E+00 


fpd(T 


1.38687E+01 


6.50911E+01 


6.17796E-01 


1.48098E+00 


3.07154E-01 


f pdcj 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+OO 


1.78179E-01 


O.OOOOOE+00 


9pd(T 


8.79879E-01 


1.03750E+00 


7.89515E-01 


8.98005E-01 


6.91073E-01 


^pdir 


-1.71792E-01 


1.25909E+00 


2.23528E+01 


1.12385E+02 


1.32243E+02 


fpd-iT 


-3.17275E-01 


-5.27343E-01 


-5.63614E+00 


-3.39863E+01 


-4.33354E+01 


f pdiT 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+OO 


6.13082E-01 


O.OOOOOE+00 


QpdTT 


7.26150E-01 


8.11349E-01 


1.09404E+00 


1.15426E+00 


1.23322E+00 


eddcT 


-3.55884E+00 


-3.25060E+01 


-1.97035E-01 


-5.40173E+02 


-2.67754E+03 


fdd(T 


8.69101E-01 


8.16200E+00 


3.64506E-01 


1.59300E+02 


6.84560E+02 


f dd(T 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+OO 


-5.69928E+00 


O.OOOOOE+OO 


gdda 


7.21756E-01 


9.54063E-01 


8.52778E-01 


1.27106E+00 


1.43875E+00 




1 6771QE+01 


-1.91191E+00 


4.03560E-01 


—2 08099E+01 


— 9.32516E— 01 


fdd-iT 


-4.08873E+00 


9.91203E-02 


-1.57338E+00 


1.75627E+00 


1.22212E-01 


f dd-n 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+OO 


-8.08524E-03 


O.OOOOOE+00 


gddn 


8.72084E-01 


7.56434E-01 


9.99799E-01 


1.10078E+00 


7.74043E-01 


BddS 


-5.68687E+01 


3.04918E+01 


3.44139E-01 


-1.01318E+02 


-4.20024E+01 


fddS 


1.99517E+01 


-3.74711E+00 


-2.87068E-02 


1.94265E+01 


8.89666E+00 


fddS 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+OO 


7.07164E-01 


O.OOOOOE+00 


gddS 


1.12436E+00 


1.06028E+00 


7.01862E-01 


1.19473E+00 


1.14476E+00 
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TABLE VIII. Tight-binding parameters (p|-|lO[) for Manganese, Iron, Cobalt, Nickel, and Copper. The equation numbers 
refer to the parametrization forms discussed in the text. 



Parameter 


Manganese 


Iron 


Cobalt 


Nickel 


Copper 


On-Sito Parameters Eqs. ( 


/) and (3) 


A 


1.65076E+00 


1.61592E+00 


1.43580E+00 


1.43847E+00 


1.47292E+00 


as 


2.82713E-01 


3.13939E-01 


3.11786E-01 


2.26020E-01 


2.13398E-02 


bs 


-2.44625E+01 


-1.50341E+02 


3.03424E+01 


2.85451E+01 


6.93939E+01 


Ca 


8.90173E+04 


5.6Q132E+Q4 


5.76787E+02 


-1.40381E+02 


-2.44787E+03 


ds 


O.OOOOOE+00 


O.QQQQQE+QQ 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+OO 


ap 


7.92094E-01 


7.50228E-01 


4.87203E-01 


4.69576E-01 


5.55994E-01 


bp 


7.18208E+01 


1.32173E+02 


4.91535E+01 


5.15499E+01 


4.82152E+01 


Cp 


3.50298E+02 


-4.83053E+03 


-8.82385E+02 


-1.16395E+03 


-3.82848E+02 


dp 


O.OOOOOE+00 


O.OOOOOE+QQ 


O.OOOOOE+00 


O.OOOOOE+OO 


O.OOOOOE+OO 


ad 


9.26905E-02 


6.79136E-Q2 


4.44373E-02 


3.58969E-02 


1.99140E-02 


bd 


1.53106E+01 


1.8Q488E+01 


1.28650E+00 


9.52349E-01 


3.46551E-01 


Cd 


9.82070E+02 


-1.71454E+03 


1.44172E+01 


-1.45504E+01 


8.08263E+01 


dd 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+OO 


O.OOOOOE+OO 






Hamiltonian Parameters 


Eq. ([lO) 






3.42882E-01 


4.55044E-01 


-1.22054E+00 


-1.38151E+00 


-9.10328E+00 


fssa 


-2.15584E-01 


-1.21734E+00 


-1.03624E-01 


-1.20825E-01 


5.17873E-01 


Qssa 


7.00302E-01 


9.14315E-01 


8.35761E-01 


8.33833E-01 


9.46662E-01 


Gspcr 


5.68705E-01 


8.72706E-01 


1.99538E+00 


2.08749E+00 


3.10063E+00 


fspcT 


1.27669E-01 


-2.65086E-03 


5.04173E-02 


8.44512E-02 


-1.14987E-01 


Qspa 


6.66959E-01 


6.78102E-01 


8.17059E-01 


8.07510E-01 


8.16824E-01 


fippa 


1.22039E+00 


1.02661E+00 


8.07087E-01 


6.39053E-01 


1.08944E+00 


fppcr 


-3.56575E-02 


4.62878E-02 


8.08439E-02 


8.01499E-02 


8.81043E-02 


Qppa 


6.20215E-01 


6.49770E-01 


6.77788E-01 


6.63093E-01 


7.44048E-01 


(^ppTT 


-9.85276E+01 


-3.61139E+01 


-1.40024E-01 


3.38998E-01 


-3.47451E+00 


fppiT 


3.70519E+01 


8.18234E+00 


-3.67928E-01 


-3.32640E-01 


4.87501E-01 


Qpp-n 


1.25798E+00 


1.17142E+00 


8.87256E-01 


8.46109E-01 


8.54245E-01 


^sda 


1.08489E+01 


5.06145E-01 


-1.26092E+00 


-1.29152E+00 


-2.32022E+00 


fsda 


-2.59008E+00 


-3.03407E-01 


8.27292E-02 


7.48429E-02 


-3.15366E-02 


gsdcr 


9.25004E-01 


8.23582E-01 


8.35781E-01 


8.30637E-01 


9.25016E-01 




5.60769E+00 


3.62443E+00 


-3.89270E-01 


-4.18197E-01 


-8.93442E-01 


fpdfT 


-1.52163E+00 


-1.20097E+00 


2.21672E-03 


-8.47911E-03 


6.29553E-02 


Qpda 


8.50729E-01 


8.70892E-01 


7.15993E-01 


6.98393E-01 


7.37891E-01 


^pdTT 


-2.24611E+00 


-1.23170E+00 


3.93950E-01 


9.56716E-02 


-5.40577E+00 


fpdn 


7.81059E-01 


8.46862E-01 


8.86907E-01 


8.22286E-01 


2.97343E+00 


Qpdn 


9.39607E-01 


9.92416E-01 


1.02442E+00 


1.03137E+00 


1.11838E+00 


eddcT 


-1.42480E+00 


-1.30200E+00 


-1.17737E+00 


-1.14634E+00 


-1.76354E+00 


fddfT 


6.29484E-03 


9.13566E-02 


-9.90018E-02 


-7.89742E-02 


1.19111E-01 


gddtj 


8.46043E-01 


7.96761E-01 


8.93869E-01 


9.02008E-01 


8.86270E-01 


^ddTT 


5.87164E+00 


3.03158E+00 


2.12933E+00 


3.62103E+00 


4.33011E+00 


fddn 


6.99707E-01 


-2.29520E-01 


6.30960E-01 


4.97609E-01 


1.60096E-01 


gddn 


1.07880E+00 


9.29910E-01 


1.04447E+00 


1.07897E+00 


1.08524E+00 


EddS 


-5.44863E+00 


-2.42866E+00 


-5.50043E-01 


-1.36239E+00 


8.61842E+00 


fddS 


-2.54667E+01 


3.43810E-01 


5.72799E-02 


2.05114E-01 


-3.97339E+00 


QddS 


1.45341E+00 


1.01267E+00 


8.74734E-01 


9.76017E-01 


1.29258E+00 
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Overlap Parameters Eq. (10) 
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7 fin4^9F— m 


O .4:^0 i Ui-J U J- 


Q OOSfi^F— 01 


Q ^1 04SF— 01 




Gppa 




7 7QnzL7F-l-nn 


9 OS^fiOF 01 


Q QQf^l /IF 01 


8 fi71 zlfiP ni 


f 

Jppa 


Zr . O U '4: _L U -L/ u u 




— ^ 7Q1 QfiF_09 


— fi 0407nF — 09 


9 n8'in4P — 03 


Qppo- 






^ S^'^OQF— 01 


^ fil '^1 9F— 01 


fi ofisiinE— 01 


^ppir 


—4 63276E+00 


—5 28179E+00 


3.22718E— 01 


4.28633E— 01 


4.59109E— 02 


fppn 


2.51612E+00 


1.65817E+00 


-8.25639E-03 


2.84019E-03 


5.34381E-02 


Qppn 


9.21078E-01 


8.02822E-01 


6.14330E-01 


5.81584E-01 


6.65480E-01 


^sda 


-1.25049E+03 


-6.29774E+02 


1.45688E+03 


1.45970E+03 


6.33227E+00 


fsdcr 


2.77503E+02 


1.71529E+02 


-3.47915E+02 


-3.44722E+02 


1.08483E+00 


Qsda 


1.41062E+00 


1.30437E+00 


1.29390E+00 


1.35652E+00 


1.12244E+00 


^pda 


-5.01714E-01 


-3.79545E+00 


-3.42786E+00 


-3.98726E+00 


-3.61230E+00 


fpda 


-9.39420E-02 


7.64597E-01 


8.06002E-01 


7.54214E-01 


2.03040E-01 


9pd(7 


7.58912E-01 


8.45424E-01 


7.56506E-01 


7.69220E-01 


9.56356E-01 




1.30596E+02 


1.21879E+02 


-1.25831E+03 


-1.47870E+03 


-3.42670E+02 


fpdn 


-3.89023E+01 


-4.37850E+01 


2.66826E+02 


2.02738E+02 


-5.22754E+01 


Qpd-K 


1.16633E+00 


1.33895E+00 


1.44349E+00 


1.46186E+00 


1.51831E+00 


&dda 


-1.50368E+03 


-4.33297E+00 


4.83403E+02 


3.09179E+03 


-1.72506E+03 


fddtj 


4.28271E+02 


4.93980E+00 


8.16729E+01 


-3.49314E+02 


4.47320E+02 


Qdda 


1.43291E+00 


1.19874E+00 


1.53933E+00 


1.55682E+00 


1.53713E+00 


GddTv 


-1.18848E+00 


-1.03641E+00 


-1.07889E+00 


-8.68984E-01 


-6.49993E-01 


fddiT 


2.24642E-02 


1.93977E-01 


1.24411E-02 


6.12811E-02 


-3.47252E-01 


Qddir 


9.42798E-01 


8.98115E-01 


8.71866E-01 


9.09043E-01 


1.08398E+00 


e-dds 


-6.61638E+01 


4.62136E+00 


-1.29889E+06 


-4.61448E+05 


2.17112E+01 


fddS 


1.62728E+01 


-1.26354E+00 


3.38458E+05 


3.24162E+05 


-4.38356E+00 


QddS 


1.15886E+00 


1.05392E+00 


1.97301E+00 


2.01457E+00 


1.18724E+00 
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TABLE IX. Tight-binding parameters (JQ^ for Strontium, Yttrium, Zirconium, Niobium, and Molybdenum. The equation 
numbers refer to the parametrization forms discussed in the text. 



Parameter 


Strontium 


Yttrium 


Zirconium 


Niobium 


Molybdenum 


On-Site Parameters Eqs. ( 


/) and (3) 


A 


1.14372E-I-00 


1.27304E+00 


9.92096E-01 


1.34109E+00 


1.35071E+00 


as 


1.77264E-01 


7.28123E-02 


4.56862E-01 


2.50727E-02 


8.66826E-02 


bs 


-1.31630E+01 


9.32506E+00 


-1.73980E+00 


5.47470E+01 


6.34195E+01 


Cs 


1.37270E+03 


-4.92331E+Q2 


-2.37835E+01 


5.11141E+03 


9.36509E+02 


ds 


O.OOOOOE+OO 


O.QQQQQE+QQ 


2.45157E+02 


O.OOOOOE+OO 


O.OOOOOE+OO 


ap 


1.76186E-01 


4.14728E-01 


3.26016E-01 


5.88564E-01 


7.20981E-01 


bp 


3.05689E+01 


7.32116E+01 


3.35603E+00 


3.39986E+01 


2.78880E+01 


Cp 


-5.40490E+02 


-1.49339E+03 


9.35236E+00 


1.02038E+03 


1.72325E+03 


dp 


O.OOOOOE+OO 


O.QQQOQE+00 


-3.24915E+01 


O.OOOOOE+OO 


O.OOOOOE+OO 


ad 


2.41720E-01 


2.Q7908E-01 


1.09229E-01 


1.38198E-01 


1.02179E-01 


bd 


1.87317E+00 


-5.61629E+00 


1.61113E+00 


5.94294E+00 


3.21725E+00 


Cd 


1.36064E+01 


1.49923E+03 


-6.46583E+00 


1.48747E+02 


1.61074E+02 


dd 


O.OOOOOE+OO 


O.OOOOOE+OO 


1.91259E+01 


O.OOOOOE+OO 


O.OOOOOE+OO 






Hamiltonian Parameters 


Eq. ([lO) 






1.01278E+02 


7.72026E+00 


5.69335E+00 


-2.50992E+01 


-7.59134E-01 


fssa 


-1.60644E+01 


-1.50419E+00 


-1.57120E+00 


-4.82208E-01 


-2.42755E+00 


Qssa 


9.17082E-01 


7.55938E-01 


8.41753E-01 


1.02758E+00 


9.86587E-01 


Gspcr 


-7.51014E+00 


-1.39670E+00 


-1.05737E+00 


1.53790E+00 


1.28784E+00 


fspcT 


1.26677E+00 


3.59513E-01 


3.53420E-01 


9.57857E-02 


2.35356E-01 


Qspcr 


7.20485E-01 


6.77441E-01 


7.03266E-01 


7.66846E-01 


7.88792E-01 


€ppcr 


-5.84498E-01 


1.02671E+00 


1.38018E+00 


1.83668E+00 


1.79716E+00 


fppcr 


1.04771E-01 


-2.50961E-02 


-8.27718E-02 


-8.85229E-02 


2.95321E-02 


Qppa 


5.41974E-01 


5.44555E-01 


5.42648E-01 


6.16705E-01 


6.71700E-01 


(^ppTT 


3.78163E+07 


9.35788E+09 


4.92810E+00 


-5.05098E+01 


-2.22510E+01 


fppiT 


-8.39757E+06 


-1.43128E+09 


1.43919E+01 


6.76702E+00 


1.82241E+00 


Qpp-n 


1.72686E+00 


2.03112E+00 


1.27121E+00 


1.15006E+00 


1.05311E+00 


^sda 


1.38973E+00 


7.11340E-01 


3.28602E-01 


-8.21822E-01 


-4.99283E-01 


fsda 


-3.46441E-01 


-1.55146E-01 


-1.96275E-01 


-6.54205E-03 


3.80005E-02 


gsdcr 


7.25100E-01 


6.31695E-01 


6.91552E-01 


7.23214E-01 


6.83391E-01 




7.31763E-01 


7.96314E+00 


5.43134E+00 


-2.56118E-01 


7.47595E-02 


fpda 


-9.56001E-02 


-1.33491E+00 


-1.54334E+00 


-1.54813E-02 


-5.07410E-02 


Qpda 


5.59815E-01 


7.43871E-01 


8.11931E-01 


6.02644E-01 


6.03739E-01 


^pdTT 


3.58801E-01 


7.43705E-02 


-3.09334E-03 


2.13205E-01 


5.98715E-01 


fpdn 


-3.27229E-02 


3.88827E-03 


2.47270E-03 


-9.42484E-03 


-2.22008E-02 


Qpdn 


5.12420E-01 


6.34856E-01 


5.14096E-01 


6.92239E-01 


7.39073E-01 


eddcT 


5.65896E+00 


-1.28617E-02 


-9.18974E-01 


-1.53648E+00 


-1.33254E+00 


fddfT 


-9.15043E~01 


-4.00283E-02 


-4.92553E-02 


2.71430E-02 


-8.68652E-03 


gddtj 


7.63696E-01 


6.20401E-01 


7.03930E-01 


7.33958E-01 


7.61081E-01 


^ddTT 


-3.38278E+00 


9.61926E-01 


1.05904E+00 


2.78488E+00 


3.48251E+00 


fddn 


9.46106E-01 


1.89330E-01 


1.63148E-01 


5.06869E-02 


-1.22045E-01 


gddn 


8.24676E-01 


7.83094E-01 


8.00179E-01 


8.57695E-01 


8.56524E-01 


EddS 


2.15625E+02 


6.79257E+01 


3.58548E+01 


-3.45690E+01 


-7.40531E-01 


fddS 


-3.46331E+01 


"1.31752E+01 


-8.49137E+00 


-1.27137E+01 


7.15057E-02 


QddS 


1.07240E+00 


1.04788E+00 


1.10660E+00 


1.29063E+00 


7.76694E-01 
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Overlap Parameters Eq. (10) 
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1 . luyoo-Cjn^uu 


f 

Jppa 




'4:. C''4:U i OJ-Jn^UU 








Qppo- 




8 01 'i43E— 01 






fi 4fifi9'iF— 01 


^ppir 


— 1 69469E+02 


—8 46463E+10 


2.23043E+01 


— 6.84291E— 01 


— 1 58174E+00 


fppn 


2.77774E+01 


1.43007E+10 


1.29974E+01 


2.91908E-01 


3.04289E-01 


Qppn 


8.77067E-01 


2.05732E+00 


1.20161E+00 


2.75732E+00 


6.45588E-01 


^sda 


-6.09261E+02 


-1.90620E+02 


-8.21444E+02 


2.67689E+00 


4.31740E-01 


fsdcr 


9.90033E+01 


3.84507E+01 


1.91013E+02 


2.37409E-01 


1.31467E-01 


Qsda 


9.57395E-01 


9.65863E-01 


1.14321E+00 


8.09336E-01 


6.84817E-01 


^pda 


-1.65574E+02 


1.29405E+00 


O.OOOOOE+00 


5.73767E-01 


4.03350E-01 


fpda 


2.72554E+01 


-9.61853E-02 


O.OOOOOE+00 


-3.90244E-02 


-3.04020E-02 


9pd(7 


8.68471E-01 


5.37486E-01 


l.OOOOOE+00 


5.08340E-01 


3.85790E-01 




6.84374E+00 


6.14466E+01 


3.66168E+01 


6.15636E-01 


4.48953E-01 


fpdiT 


-1.43215E+00 


-1.25963E+01 


-9.03817E+00 


-2.16620E-01 


-1.59315E-01 


QpdTV 


7.16925E-01 


9.14935E-01 


9.83872E-01 


7.64261E-01 


7.32118E-01 


&dda 


-1.39363E+01 


-3.24712E+01 


-4.84660E+01 


-3.43653E-01 


2.47289E-01 


fddtj 


2.65537E+00 


6.51211E+00 


1.17934E+01 


1.55909E-01 


-1.72562E-02 


Qdda 


7.39857E-01 


8.29440E-01 


1.02194E+00 


6.65039E-01 


4.36793E-01 


Gdd-K 


1.37850E+02 


1.44898E+02 


1.56704E+02 


-1.88524E+00 


-7.55891E-01 


fddn 


-2.52973E+01 


-2.91678E+01 


-3.85906E+01 


2.56521E-02 


5.55713E-02 


9dd7T 


9.08764E-01 


9.71214E-01 


1.10151E+00 


7.54927E-01 


6.03483E-01 


GddS 


-3.88013E+02 


8.80037E+01 


-1.14800E+01 


2.48974E-01 


1.21646E+00 


fddS 


8.01182E+01 


4.30951E+00 


5.38231E+00 


-2.44331E-02 


6.30985E-03 


QddS 


1.10480E+00 


1.15554E+00 


1.11935E+00 


5.59223E-01 


8.82482E-01 



17 



TABLE X. Tight-binding parameters (p|-|lC|) for Technetium, Ruthenium, Rhodium, Palladium, and Silver. The equation 
numbers refer to the parametrization forms discussed in the text. 



Parameter 


Technetium 


Ruthenium 


Rhodium 


Palladium 


Silver 


On-Site Parameters Eqs. ( 


1) and (3) 


A 


1.35605E+00 


1.34253E+00 


1.34056E+00 


1.27733E+00 


1.22495E+00 


as 


1.56969E-01 


8.93160E-02 


9.12651E-02 


9.21895E-02 


1.71377E-01 


bs 


4.96528E+01 


4.12031E+01 


4.48975E+01 


2.66437E+01 


8.37704E+00 


Cs 


6.87483E+02 


1.29943E+Q3 


8.96404E+02 


2.45427E+02 


-1.77374E+01 


ds 


O.OOOOOE+OO 


O.QOOQQE+QQ 


O.OOOOOE+OO 


O.OOOOOE+OO 


O.OOOOOE+OO 


ap 


7.27706E-01 


6.70655E-01 


6.20365E-01 


5.05324E-01 


6.05047E-01 


bp 


2.33596E+01 


4.08917E+01 


4.44309E+01 


3.52846E+01 


1.23316E+01 


Cp 


1.80513E+03 


3.65045E+01 


-7.16816E+01 


-1.73870E+02 


2.00392E+01 


dp 


O.OOOOOE+OO 


O.QQQOQE+QQ 


O.OOOOOE+OO 


O.OOOOOE+OO 


O.OOOOOE+OO 


ad 


8.87546E-02 


6.47777E-02 


4.90478E-02 


2.77484E-02 


4.32432E-03 


bd 


1.68547E+00 


1.15434E+00 


7.99566E-01 


2.34504E-01 


8.54569E-02 


Cd 


9.14454E+01 


4.02940E+01 


-1.33019E+01 


-1.16871E+01 


-5.98728E-01 


dd 


O.OOOOOE+OO 


O.OOOOOE+OO 


O.OOOOOE+OO 


O.OOOOOE+OO 


O.OOOOOE+OO 






Hamiltonian Parameters 


Eq. ([lO) 






-2.42325E-01 


-1.06826E+01 


-1.03870E+01 


-9.65803E+00 


-3.02245E+00 


fssa 


-2.35701E+00 


-5.19162E-01 


O.OOOOOE+OO 


-6.58566E-02 


2.82192E-01 


Qssa 


9.90587E-01 


9.98652E-01 


9.74343E-01 


9.65567E-01 


7.53400E-01 


Gspcr 


1.32706E+00 


5.36202E+00 


7.11658E+00 


6.92889E+00 


2.46216E+00 


fspcT 


2.33454E-01 


2.42663E-02 


O.OOOOOE+OO 


-8.05132E-03 


-6.72462E-02 


Qspcr 


7.91274E-01 


8.79989E-01 


9.06080E-01 


8.93397E-01 


7.93558E--01 


Gppfj 


1.67313E+00 


3.00393E+00 


2.61732E+00 


3.10178E+00 


-1.52677E+01 


fppcr 


7.40188E-03 


1.98109E--02 


O.OOOOOE+OO 


-3.64082E-03 


4.55095E+00 


Qppa 


6.65741E-01 


7.26256E-01 


7.04652E-01 


7.29562E-01 


8.98600E-01 


(^ppTT 


-2.05304E+01 


-8.52581E-02 


-9.15155E-02 


-8.63018E-02 


-4.24238E+02 


fppiT 


2.07455E+00 


-2.22726E-03 


O.OOOOOE+OO 


-3.08492E-03 


-2.33808E+02 


Qpp-n 


1.03629E+00 


5.41216E-01 


5.24756E-01 


5.72566E-01 


1.61153E+00 


^sda 


-6.50753E-01 


-1.63451E+00 


-4.27539E+00 


-1.79028E+00 


-1.12908E+00 


fsda 


5.67702E-02 


-9.54826E-02 


O.OOOOOE+OO 


-1.40935E-02 


9.41454E-02 


gsdcr 


6.57744E-01 


8.93940E-01 


9.40417E-01 


8.55849E-01 


7.06503E-01 




-2.26321E-02 


-5.08077E-01 


-6.91759E-01 


-7.13093E-01 


3.72317E+00 


fpdfT 


-5.91240E-02 


6.17604E-04 


O.OOOOOE+OO 


-2.66134E-03 


-1.48971E+00 


Qpda 


6.24648E-01 


6.66661E-01 


6.72844E-01 


6.91813E-01 


9.00056E-01 


^pdTT 


5.73187E-01 


3.33146E-01 


1.05800E-01 


2.53192E-01 


2.76295E+00 


fpdn 


-3.81543E-02 


-3.86754E-03 


O.OOOOOE+OO 


-2.89608E-03 


1.87799E~01 


Qpdn 


7.46410E-01 


7.10996E-01 


6.16605E-01 


7.19375E-01 


1.02993E+00 


eddcT 


-1.32507E+00 


-2.47448E+00 


-2.24220E+00 


-2.00840E+00 


-2.83314E+00 


fddfT 


-9.81573E-03 


-4.63223E-02 


O.OOOOOE+OO 


-1.50392E-02 


-6.41209E-01 


gdda 


7.70176E-01 


8.54974E-01 


8.45401E-01 


8.53056E-01 


9.78615E-01 


^dd-K 


4.14667E+00 


4.17666E+00 


5.70876E+00 


8.17277E+00 


9.65404E+00 


fddn 


-1.30743E-01 


1.94391E-02 


O.OOOOOE+OO 


-8.64748E-02 


-8.32773E-01 


QddiT 


8.91789E-01 


9.34693E-01 


9.82160E-01 


1.02752E+00 


1.01612E+00 


EddS 


-9.72649E-01 


5.28412E+00 


-1.68497E+00 


-2.34604E+00 


2.59720E+01 


fddS 


1.30476E-01 


~2.14644E+00 


O.OOOOOE+OO 


3.25537E-01 


-1.38837E+01 


QddS 


7.84398E-01 


1.11033E+00 


1.04917E+00 


9.76385E-01 


1.36049E+00 
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Overlap Parameters Eq. (10) 
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1 40686E+01 

_L • j: W W Vj W J — J 1 W -L 


— 1.40201E— 02 


— 6.16862E— 01 


—5 78210E+03 


fppn 


2.88300E-01 


-9.07707E-04 


O.OOOOOE+00 


1.33749E-01 


1.25726E+03 


Qpp-K 


6.50767E-01 


1.98112E+00 


1.73125E+00 


7.11976E-01 


1.36538E+00 


GsdcT 


7.28149E-01 


2.66680E-01 


-2.09986E+00 


-1.37473E+01 


2.87589E-01 


fsdcr 


7.43369E-02 


1.19878E-01 


O.OOOOOE+00 


2.65050E-01 


3.25719E-02 


Qsda 


6.87723E-01 


7.69449E-01 


1.16939E+00 


1.07480E+00 


7.24168E-01 


&pda 


4.34355E-01 


-1.54734E+03 


-4.40618E+02 


-1.25579E+02 


-1.24324E+01 


fpdo- 


-3.40641E-02 


-2.04299E+03 


O.OOOOOE+00 


-3.20496E+01 


1.90544E+00 


Qpda 


4.18094E-01 


1.95694E+00 


1.38332E+00 


1.25500E+00 


9.52608E-01 


GpdTT 


5.61187E-01 


-3.96833E-01 


-5.68639E+01 


-7.57995E+02 


-1.79373E+00 


fpdiT 


-1.39188E-01 


5.43342E-02 


O.OOOOOE+00 


-3.66651E+00 


2.67272E-01 


Qpd-K 


7.25597E-01 


5.68958E-01 


1.21849E+00 


1.37921E+00 


7.56593E-01 


Cdda 


3.23273E-01 


-4.42200E+00 


2.38357E+00 


3.35457E+01 


1.43164E+01 


fddtj 


1.72275E-02 


2.03106E+00 


O.OOOOOE+00 


-6.81529E-01 


-2.47557E+00 


Qdda 


6.14412E-01 


9.65997E-01 


8.81490E-01 


1.14453E+00 


9.44667E--01 


GddTT 


-6.98334E-01 


-4.84435E+00 


-2.32596E+01 


-1.50866E+01 


-7.84762E+00 


fddn 


5.95975E-02 


-4.04300E-01 


O.OOOOOE+00 


1.44912E-01 


2.94062E-01 


9dd7T 


6.29126E-01 


1.02395E+00 


1.12774E+00 


1.10093E+00 


1.04716E+00 


CddS 


1.39075E+00 


1.33535E+04 


5.83229E+03 


1.05067E+03 


5.03416E+01 


fddS 


-5.99108E-02 


4.97401E+03 


O.OOOOOE+00 


2.44615E+01 


-7.80465E+00 


gddS 


9.31751E-01 


1.76428E+00 


1.60401E+00 


1.47110E+00 


1.15950E+00 
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TABLE XI. Tight-binding parameters for Barium, 

to the parametrization forms discussed in the text. 



Hafnium, Tantalum, and Tungsten. The equation numbers refer 



Parameter 


Barium 


Hafnium 


Tantalum 


Tungsten 


On-Site Parameters Eqs. ( 


() and ( 


9) 



A 


1.11634E+00 


1.19630E+00 


1.27675E+00 


1.32674E+00 


as 


1.50138E-01 


9.10131E-01 


3.19292E-01 


4.33397E-01 


bs 


-1.53745E+01 


5.29294E+00 


3.66642E-01 


1.42840E+00 


Ca 


1.56954E+03 


5.04222E+02 


4.42664E+02 


4.00732E+02 


ds 


O.OOOOOE+00 


4.95952E+02 


O.OOOOOE+00 


O.OOOOOE+00 


ap 


1.28894E-01 


4.65965E-01 


3.59226E-01 


4.27055E-01 


bp 


2.72310E+01 


2.94282E+01 


5.38397E+01 


8.21465E+01 


Cp 


1.64446E+03 


-3.86257E+02 


-4.57578E+02 


-1.13769E+03 


dp 


O.OOOOOE+00 


1.23293E+03 


O.OOOOOE+00 


O.OOOOOE+00 


ad 


1.82919E-01 


2.63519E-01 


1.61383E-01 


2.17265E-01 


bd 


9.59617E-01 


4.75514E+00 


2.74648E+00 


-2.63465E+00 


Cd 


-9.73190E+01 


4.74879E+01 


1.22440E+02 


3.23869E+02 


dd 


O.OOOOOE+00 


-7.45384E+02 


O.OOOOOE+00 


O.OOOOOE+00 


Hamiltonian Parameters Eq. (10) 




1.07625E+02 


9.13877E+00 


-3.81526E-01 


4.59797E+00 


fssa 


-1.54756E+01 


1.54875E+00 


O.OOOOOE+00 


-1.06279E+00 


Qssa 


8.8354 7E-01 


9.83725E-01 


5.91644E~01 


7.81632E-01 


Gspa 


-8.15023E+00 


-9.83656E-01 


1.07250E+00 


-2.71873E-02 


fspcT 


1.18680E+00 


4.73735E-01 


O.OOOOOE+00 


4.18002E-01 


Qspa 


7.05573E-01 


6.25311E~01 


6.60713E-01 


7.28960E-01 




-7.65354E-01 


-3.62500E+00 


1.96224E+00 


1.04741E+00 


fppcr 


1.11812E-01 


1.28779E+00 


O.OOOOOE+00 


1.14611E+00 


Qppa 


5.22978E-01 


7.57257E-01 


6.54733E-01 


7.99230E-01 


f^ppTT 


-4.49769E+08 


-1.86311E+00 


-1.25096E+02 


4.35429E+01 


fppTT 


6.24773E+07 


1.55799E-01 


O.OOOOOE+00 


-7.09146E+00 


QppiT 


1.72686E+00 


9.06562E-01 


1.33474E+00 


1.10993E+00 


^sda 


1.51803E+00 


3.28266E+00 


-5.39923E-01 


-7.48675E-01 


fsda 


-3.17484E-01 


-1.20838E+00 


O.OOOOOE+00 


-5.87424E-02 


QadtT 


7.01455E-01 


7.79418E--01 


7.12949E-01 


7.15275E-01 


^pda 


1.38944E+00 


1.08653E-01 


-1.52215E-01 


5.22978E+00 


fpdfT 


-1.70619E-01 


-5.64207E-01 


O.OOOOOE+00 


-1.54532E+00 


Qpda 


5.57474E-01 


7.13632E-01 


5.71662E-01 


8.11224E-01 


^pdTT 


3.19179E-01 


-3.46529E+00 


3.74952E-01 


-3.92986E+01 


fpd-R 


-3.37340E-02 


9.67868E-01 


O.OOOOOE+00 


7.65413E+00 


Qpdn 


6.20946E-01 


1.05238E+00 


6.53815E~01 


1.03638E+00 


eddcr 


9.48343E+00 


-2.15888E-01 


-1.30670E+00 


-6.77733E-01 


fdda 


-1.44518E+00 


-2.80122E-01 


O.OOOOOE+00 


-5.81468E-01 


Qdda 


7.49400E-01 


7.11762E-01 


7.25180E-01 


8.48777E-01 


^ddTT 


-2.85791E+00 


1.71474E+00 


1.12618E+00 


5.95483E+00 


fdd-R 


8.93926E-01 


5.62278E-01 


O.OOOOOE+00 


-3.94918E-01 


gddn 


8.11059E-01 


9.14855E-01 


7.22741E-01 


8.85080E-01 


eddS 


1.87063E+02 


1.33656E+01 


-2.46591E+00 


-3.0184 7E+01 


fddS 


-2.53043E+01 


6.13820E-01 


O.OOOOOE+00 


2.76809E+00 


gddS 


9.83611E-01 


1.07864E+00 


9.37069E-01 


1.19769E+00 
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Overlap Parameters Eq. (LO) 





1 S'^799F'-U-ni 

1 .00 / ZZr£!j^U-L 


1 189'^lF4-ni 


A Q'i9QfiF-l-nn 


1 ni 1 SfiF-KT^ 


f 

J ssa 








9 fi4fi9QF-l-n9 


Qssa 


( . _L C/ ^ ^ J—J \J i- 


Q 4fi'i83R— ni 


8 'i9fi9flF,— ni 


1 9'i3fl7F,-l-00 


Gsp<j 






A 7'^^78F4-n9 


Q 7Sc:7QF m 


f 

J spa 








— 1 OQ^fil F— m 


Qspcr 


1 MQDSR-i-nr) 






7 9fi1 3SF,— ni 


Gppa 


7 '^'^1 '^nF4^n9 

i .00 i-0\jrj~r\jzi 


1 '^4n'^i F-i~ni 

-L .^J'y:UO-L-L!J^^U-L 




7 Q4Q'^7F-l-ni 


f 

Jppa 






n nnnnnF,-i-nn 




Qppo- 






1 9Q4S9F-l-nn 


1 1 97q'iF,-l-00 


^ppir 


—2 62750E+01 


2.80432E+01 


1 70576E+04 


2.84202E+02 


fppn 


1.77696E+01 


-2.12441E+00 


O.OOOOOE+00 


-4.98337E+01 


Qpp-K 


9.13209E-01 


1.07278E+00 


2.02320E+00 


1.15218E+00 


^sda 


-5.54362E+02 


6.14493E+00 


2.10081E+00 


1.86528E+01 


fsdcr 


8.25381E+01 


-1.61559E+00 


O.OOOOOE+00 


1.45875E+01 


Qsda 


9.13548E-01 


1.28745E+00 


7.03679E-01 


1.19625E+00 


^pda 


-1.98494E+02 


7.95537E+00 


1.08905E+00 


-2.31625E+01 


fpda 


3.13853E+01 


-2.35268E+00 


O.OOOOOE+00 


1.40867E+01 


9pd(7 


8.78901E-01 


9.07766E-01 


6.11045E-01 


1.14689E+00 




9.42421E+00 


1.55087E+00 


-3.47807E-02 


-2.20429E+01 


fpdir 


-1.77293E+00 


-4.88415E-01 


O.OOOOOE+00 


1.06219E+00 


Qpd-K 


6.93194E-01 


8.56422E-01 


4.74891E-01 


1.06161E+00 


Cdda 


-2.28081E+01 


3.36607E+01 


2.24207E-01 


-4.53347E+01 


fddtj 


3.96599E+00 


-8.74169E+00 


O.OOOOOE+00 


1.19568E+01 


Qdda 


7.50292E-01 


1.20766E+00 


4.92531E-01 


1.03146E+00 


GddTT 


1.30794E+02 


-1.26460E+01 


-1.83289E+00 


5.24024E+02 


fddn 


-2.48946E+01 


1.05554E+00 


O.OOOOOE+00 


-1.33112E+02 


Qddn 


8.94081E-01 


9.07237E-01 


7.62263E-01 


1.22976E+00 


edds 


1.07020E+04 


6.04877E+01 


4.42902E-01 


-1.38362E+01 


fddS 


-1.18878E+03 


-5.00466E+00 


O.OOOOOE+00 


6.70560E+00 


QddS 


1.17919E+00 


1.09546E+00 


7.43396E-01 


1.19251E+00 



21 



TABLE XII. Tight-binding parameters ([7|-|lC|) for Rhenium, Osmium, Iridium, Platinum, and Gold. The equation numbers 
refer to the parametrization forms discussed in the text. 



Parameter 


Rhenium 


Osmium 


Iridium 


Platinum 


Gold 


On-Site Parameters Eqs. ( 


1) and (3) 


A 


1.41690E+00 


1.43787E+00 


1.48862E+00 


1.48637E+00 


1.37357E+00 


as 


1.25774E-01 


2.19235E-01 


2.58676E-01 


5.83156E-03 


9.52842E~02 


bs 


4.32410E+01 


2.32773E+01 


8.99491E+01 


2.05116E+02 


4.03569E+01 


Cs 


-2.65528E+03 


-1.49716E+02 


-5.38470E+03 


-2.04105E+04 


-3.05963E+03 


ds 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+00 


ap 


6.08286E-01 


7.15539E-01 


6.49919E-01 


8.17823E-01 


5.75191E-01 


bp 


5.39672E+01 


7.50020E+01 


1.26272E+02 


5.89750E+01 


1.40526E+01 


Cp 


5.88296E+03 


4.09052E+03 


9.65187E+03 


1.12390E+04 


-8.66127E+00 


dp 


O.OOOOOE+OO 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+00 


ad 


1.15150E-01 


9.48499E-02 


5.69648E-02 


4.98288E-02 


2.16623E-02 


bd 


3.02548E-I-00 


-8.99667E-01 


2.78655E-01 


-4.47348E+00 


1.02921E+00 


Cd 


5.66609E+02 


5.91858E+02 


1.80727E+02 


6.15745E+02 


-3.77491E+01 


dd 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+00 


O.OOOOOE+00 






Hamiltonian Parameters 


Eq. ([lO) 






1.53798E+01 


9.56167E+00 


-1.25321E+00 


-1.38285E+00 


-3.80906E+00 


fssa 


-3.57908E+00 


-3.86825E+00 


-1.06638E-01 


-1.20674E-01 


-4.38467E-01 


Qssa 


8.70774E-01 


9.68851E-01 


8.16908E-01 


8.13945E-01 


9.08853E-01 


Gspcr 


9.12123E-01 


7.69439E-01 


1.84853E+00 


1.83354E+00 


1.15267E+00 


fspcT 


2.73036E-01 


2.88606E-01 


5.84833E-01 


5.86127E-01 


-1.43508E-02 


Qspcr 


7.96353E-01 


7.87669E-01 


8.45321E-01 


8.45134E-01 


6.80672E-01 


Gppfj 


2.22720E+00 


1.60364E+00 


1.10831E+00 


1.26085E+00 


-3.33685E+00 


fppcr 


2.72264E-01 


4.70495E-01 


1.03256E+00 


1.03890E+00 


1.05126E+00 


Qppa 


8.00001E-01 


7.56961E-01 


8.31571E-01 


8.31225E-01 


7.97642E-01 


^ppTT 


-1.26133E+01 


9.34510E-01 


1.91089E+02 


2.05176E+02 


4.78096E+02 


fppTT 


-2.52237E~01 


-4.41535E-01 


-4.16812E+01 


-3.96380E+01 


-1.14000E+02 


QppiT 


1.00962E+00 


8.42600E-01 


1.13693E+00 


1.13165E+00 


1.39758E+00 


^sda 


-6.80815E-01 


-1.20891E+00 


-3.02142E+00 


-2.89295E+00 


-2.73600E+00 


fsdcr 


5.60338E-02 


1.00979E-01 


3.22373E-01 


2.44946E-01 


2.21303E-01 


gsdcr 


5.97382E-01 


6.87352E-01 


7.78863E-01 


7.98201E-01 


7.86845E-01 




-3.89432E-01 


2.91611E-01 


1.91293E+00 


1.58551E+00 


3.72568E+00 


fpda 


-5.46612E-02 


-1.78914E-01 


-6.34580E-01 


-6.94619E-01 


-1.48066E+00 


Qpda 


7.33599E-01 


7.31519E-01 


8.28540E-01 


8.07844E-01 


8.83740E-01 


^pdTT 


6.87010E-01 


8.49872E-01 


1.12973E+00 


8.60594E-01 


-4.98763E+00 


fpd-R 


7.97327E-02 


9.73177E-02 


-5.49139E-02 


-3.28210E-02 


1.61999E+00 


QpdTT 


8.44466E-01 


8.65611E-01 


8.13738E-01 


8.59014E-01 


9.86804E~01 


eddcT 


-2.05734E+00 


--1.84915E+00 


-1.69848E+00 


-1.75160E+00 


-3.20091E+00 


fdda 


-7.81615E-02 


-9.78854E-02 


-3.17559E-01 


-3.17277E-01 


2.80249E-01 


Qdda 


8.23143E-01 


8.19999E-01 


8.67395E-01 


8.74769E-01 


8.01174E-01 


^ddTT 


4.93020E+00 


4.81023E+00 


5.45277E+00 


7.04207E+00 


1.28200E+01 


fddn 


-8.00019E-02 


-6.04252E-02 


-3.36078E-02 


-3.61063E-01 


-1.27555E+00 


gddn 


8.98430E-01 


9.11958E-01 


9.38700E-01 


9.55738E-01 


9.96619E-01 


EddS 


-3.42830E+00 


-2.11302E+00 


-8.81866E-01 


-1.19555E+00 


1.54046E+03 


fddS 


3.08896E-01 


1.76612E-01 


1.01129E-01 


1.68460E-01 


-3.98715E+02 


gddS 


9.39165E-01 


9.26725E-01 


8.09760E-01 


8.73327E-01 


1.50519E+00 
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Overlap Parameters Eq. (10) 
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f 
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Qppo- 
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G'ppTT 
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4.66106E— 01 
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—4 68566E+02 


7 25807E— 01 


fppn 


9.07266E-02 


-2.97036E-02 


-1.37636E+02 


1.29049E+02 


2.46341E-04 


Qppn 


9.40898E-01 


1.00008E+00 


1.22085E+00 


1.24801E+00 


5.60606E-01 


^sda 


6.58787E+00 


-2.99936E-01 


-1.23356E+00 


-1.88028E+00 


7.82884E+02 


fsdcr 


5.25043E-01 


6.59351E-01 


4.96759E-01 


4.33068E-01 


-5.72404E+01 


Qsda 


9.32724E-01 


7.95665E-01 


7.36238E-01 


7.48675E-01 


1.35691E+00 


^pda 


1.51040E+00 


8.30255E-01 


8.06687E-01 


1.40632E-01 


-1.59724E+01 


fpdo- 


8.03788E-02 


-6.94731E-02 


-5.95955E-02 


-1.08624E-02 


3.65486E+00 


9pd(7 


6.91649E-01 


5.55561E-01 


5.19117E-01 


4.30680E-01 


9.11474E-01 


Gpdiv 


2.77713E+00 


-1.45265E-01 


-1.75883E+00 


-2.32135E+00 


--2.09876E+03 


fpd-K 


-5.08983E-01 


-7.56943E-03 


3.88548E-01 


3.29353E-01 


3.24572E+02 


Qpdir 


7.90651E-01 


l.OOOUE+OO 
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7.81550E-01 


1.36489E+00 


&dda 


2.08835E-01 


-2.14857E+01 


8.54708E-01 


7.85665E-01 


1.16130E+04 


fddtj 


7.98506E-02 


8.12810E+00 


-7.41190E-03 


-4.01901E-02 


-6.39800E+02 


gdda 


6.28966E-01 


1.04606E+00 


6.96686E-01 


7.02103E-01 


1.57650E+00 


Gdd-K 


-1.45230E+00 


-1.80998E+00 


-3.39400E+00 


-3.26443E+00 


-2.12382E+00 


fddiT 


1.43697E-01 


1.38350E-01 


-2.60852E-02 


5.45025E-02 


2.85113E-01 


Qddir 


6.78231E-01 


7.61277E-01 


8.93060E-01 


9.05926E-01 


7.62212E-01 


e-dds 


9.19531E-02 


1.62148E-01 


2.14814E+02 


1.65627E+03 


2.32941E+01 


fddS 


7.24211E-02 


-8.40874E-03 


1.62838E+01 


-2.95917E+02 


-4.22448E+00 


QddS 


8.78356E-01 


7.18986E-01 


1.40077E+00 


1.32876E+00 


1.01272E+00 
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